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THE DEVELOPMENT AND PRELIMINARY APPLICATION OF AN 

INVARIANT COUPLED DIFFUSION AND CHEMISTRY MODEL 

Glenn R. Hilst ■ 

Coleman duP. Donaldson 
Milton Teske 
Ross Conti llano 
Johnny Freiberg 

Aeronautical Research Associates of Princeton, Inc. 

INTRODUCTION 


During the period March 1972 to March 1973 > the Aeronautical 
Research Associates of Princeton (A.R.A.P.) undertook the develop- 
ment of invariant models (second-order closure) of diffusion and • 
chemical reactions in turbulent flow systems. Several applica- 
tions of such models were evident and as a result the development 
phase of this work was jointly funded by NASA, EPA, and, by 
transfer of funds to NASA, by DOT under Contracts NASI -11433 } 

EPA 68-02-0014, and NAS1-11873. Credit for the support of the 
basic technological development of the models reported here is 
therefore shared by these three agencies. 

The applications of the models varies considerably among the 
sponsoring agencies, however, with EPA's interest centered on 
pollution problems of the atmospheric boundary layer and DOT'S 
interests centering primarily on the behavior of the far wake of 
SST exhaust products at cruise altitudes ( ~ 20 km). The problem 
orientation of our contract with NASA, NAS1-11433, centers on the 
behavior of the SST exhaust plume during its first hour or so of 
existence. One of the objectives of this study has been to 
provide guidance for the design of a possible NASA program for 
in situ measurements of SST exhaust products at cruise altitudes 

and, eventually, to use the models for data analysis and 

interpretation. ■ , . _ . • • ... . 

At the time this work was initiated a general assessment of 
the potential importance of" inhomogeneous mixing on chemical 
reaction . rates, and the. basic approach, to modeling these effects 
via second-order closure. techniques had been developed .by an 



in-house program at A.R.A.P.* At the same time development of 
the invariant diffusion models had progressed to the point where 
it was clear that a coupled diffusion/chemistry simulation system 
depended primarily on successful closure of the chemical kinetic 
equations. This report therefore directs attention primarily to 
the chemical sub-model; the reader who wishes to familiarize him- 
self with the details of the diffusion modeling is referred to a 
recent Monograph prepared by Donaldson on this topic**. Familiar 
ity with the techniques of invariant modeling of turbulence and 
diffusion is assumed in the present report. 


SYMBOLS 


A, 


B, 


Br 


F 

g 

k 


M 


n. 

i 


N 


moment term defined by Eq. (24) 

ratio of mean concentrations U /C T 

... of al 

ratio of reaction rates 

concentration of 1th chemical species 
moment term defined by Eq . (26) 
horizontal fl.ux defined by Eq. (64) 
acceleration due to gravity 



reaction rate constant 
•value of 


C-C- — p 

when C ' QX = 0 


C C ft 
a p 


a 


frequency distribution of ith species 
2 


i n i 


Coleman duP. Donaldson and Glenn R. Hilst , "Effect of 
Inhomogeneous Mixing on Atmospheric Photochemical Reactions , 1 
Environmental Science and Technology , Vol. 6, September 19/2. 

Coleman duP. Donaldson, "On the Production of Atmospheric 
Turbulence and the Dispersal of Atmospheric Pollutants, " 
A.R.A.P. Report No. 186, Vol. 1, prepared under EPA 
Contract 68-02-0014, December 1972. 
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p 

pressure 


q 

square root of twice the turbulent kinetic energy 


r 

correlation coefficient defined by Eq. (15) 


S 

wind shear 


t 

time 



intermittency factor- 


T 

absolute temperature 


T o 

adiabatic temperature 


U 

velocity 


X 

axial coordinate 


Z 

vertical distance 


X 

micro-scale length 


A i 

length scale 


A 2 

length scale 


A 3 

length scale 


v o 

kinematic viscosity 


a 

standard deviation 


z 

Superscripts 


— 

mean component 


1 

fluctuating component 


Subscripts 

chem 

reaction rate due to chemistry 


I 

reaction rate neglecting third-order correlations 
(Table 1) 


s 

steady state value (see e.g. Eq. (12) 

* 

6 

chemical species 

3 


o 


initial value 


1,2 species designation 


THE PROBLEM 


The objective of this program, and therefore the primary 
problem which we have addressed, has been to fabricate a useful 
coupled model which can simulate the combined effects of turbulent 
diffusion and chemical depletion on the concentration patterns of 
reactive chemical species emanating from common or separate 
sources. However, against the background of assessment and 
modular model development described in the previous section, it 
was evident at the initiation of the present program that two 
major problems had to be solved first. 

•1. More rigorous analyses were required in order to 
determine the magnitude of the effects of inhomogeneous mixing 
on chemical reaction rates, the conditions under which these 
effects could be realized, and an evaluation of the likelihood 
that these conditions actually occur in atmospheric pollutant 
situations. For example, if it could be shown that these effects 
were either always insignificant or constituted only a transient 
perturbation of the chemical kinetic rates, a coupled diffusion/ 
chemistry model could be readily constructed using the convention- 
al mean-value chemical kinetic equations. 

2. Given that the results of the above analyses were not 
totally negative, i.e., negative in the sense that no important 
real-world situations could be found in which concentration 
fluctuations played a significant role, it was recognized that 
the second immediate problem was the development of a useful 
closure scheme for the third-order correlations inherent in the 
complete chemical kinetic equations. It was also recognized that 
under some circumstances the third-order correlations could be 
neglected. However, this assumption restricts strongly the range 
of joint frequency distributions of reactant concentrations which 
can be considered, and denies any semblance of generality in the 
chemical sub-model. A more appropriate, although approximate, 
closure scheme was required. 

With these considerations in mind, first efforts were devoted 
to these two problems. By mid-November 1972 both had been 
resolved, and attention was focussed on the assembly of the first 
coupled diffusion/chemistry model. The results of the earlier 
work on the analyses of the magnitude and significance of inhomo- 
geneous mixing on chemical reactions and the development of a 
closure scheme at the level of third-order correlations of 
concentration fluctuations have been assembled as a technical 
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paper which was presented at the 11th Aerospace Sciences Confer- 
ence of the AIAA in Washington, D. C., January 10, 1973*. A 
slightly modified version of this paper is included in the present 
report. The major results discussed there are: 


1. There are indeed real-world atmospheric pollution 
problems in which neglect of the fluctuations of concentrations 
of reacting species introduces significant errors. 

2. An approximate closure scheme for the chemical sub-model 
which conforms to the principles of invariant modeling and which 
accounts for the effects of inhomogeneous mixing over a wide 
range of conditions (concentration variance -to-mean-squared ratios 
up to 100) has been developed. This sub-model predicts reaction 
rates to within a factor of two of the exact chemical kinetic 
solutions for situations where the mean-value chemical kinetic 
approximation incurs errors of a factor of 100. On the other hand, 
the chemical kinetic sub-model recovers the mean-value approxima- 
tion when the concentration fluctuations are indeed insignificant 
in chemical reaction rates. This second-order closure model may 
therefore be considered as a generalized (but still approximate) 
solution of the chemical kinetics equations. 


AN EVALUATION OP THE EFFECTS 
OF INHOMOGENEOUS MIXING 


The Basic Chemical Kinetic Equations for 
Inhomogeneous Mixtures 

Following Donaldson and Hilst** we assume an isothermal, 
irreversible, two-body reaction between chemical species a and 
f3 to form y and 6 . 

a + p — ► 7 + 6 (l) 

Further, we assume that the reaction rate for any joint values of 
the concentrations of the reacting chemical species are correctly 
specified by 


Glenn R. Hilst, "Solutions of the Chemical Kinetic Equations for 
Initially Inhomogeneous Mixtures," AIAA Paper No. 73 - 101, 

January 1973, Washington, D. C. 
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a 

6t 

" k 1 C C ft 

lap 

(2) 

and 


bCp 

- k 2 c a c p 

(3) 


" It (C 7 

+ °S> = 

It (C a + V 

(4) 


where denotes the concentration of the 1th chemical species 

(expressed as a mass fraction), and k]_ and k 2 are the reaction 
rate constants. 


Equations (2) and (3) specify the local instantaneous rate of 
change of the concentration of the reactants. In order to deter- 
mine the average rate of change, we assume the history of the joint 
values of C a and Cg at a fixed location comprises a stationary 
time series and each may be dissected into its mean and fluctuating 
components 


and 


C = C + C » 

(5) 

a a a 

C P = C p + C p 

(6) 


and, by definition, = Cg = 0 . Under these assumptions the 

chemical kinetic equations for the average rates of change of the 
concentrations of a and (3 at that location are readily shown 
to be 


and 


be 

a 


a 

bt - 

- k, (C C ft + C 'Cl) 
1 v a P a (3' 

(7) 

6t' 

- k 0 (C C ft + r'cT) 
2 v a P a p' 

(8) 


In order 
equation for 


to so lve Equations (J) and (8), we need a prediction 
. This is readily derived as* 



Donaldson and Hils 


t. 


op. cit.) 
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( 9 ) 


Be 'C ' _ — p p 

■ - = - k, (C CX + C R C 'C ' + C'CA) 
at 1 v a p P a P a P ' 


,2 


- k 0 (C ft C' + C C'C' + C' C') 
2 v P a aap a P' 


2 2 2 2 

which introduces four new terms, C^ , C^ , CC^ ’ and Cf 


2 2 

The prediction equations for and are 


'a w p 


BC ' c 
a 


2 


- 2k, (C 0 C + C C'C’ + C' C') 
1 ' (3 a aap a P' 


and 


( 10 ) 


Be a 


P 


-vr — = - 2k 0 (C CA 2 + C R C'C' + C'C, 
at 2 ' a f3 Pap a f 


(id 


and they do not introduce any more new terms. In order to close 
Equations (7) through (ll), and thereby achieve a chemical sub- 
model for reactions in inhomogeneous mixtures, we require predic- 

P o 

tion equations for the third-order correlations C ' CA and C 'CA . 

a p a p 

Before proceeding to the closure problem, however, it is 
instructive to examine more closely the limits of the effects of 
concentration fluctuations on chemical reaction rates and the 
conditions under which these effects become significant. This 
examination may be made in two steps, l) when are the fluctuations 
of concentration negligible (i.e., when are the reaction rates 
predicted satisfactorily by the mean values of concentration 
alone?) and 2) when may the third -order correlations be neglected? 
For the latter cases. Equations (7) through (ll) comprise the 
closed set discussed by Donaldson and Hilst*. 


The Limits of Errors in Reaction Rate Predictions 
If Concentration Fluctuations are Neglected 

Since the neglect of concentration fluctuations in determin- 
ing reaction rates is equivalent to the assumption that the local 
values of C a and Cp are constant in time, C^ = Cg = 0 and 
the reaction rates predicted under this assumption are simply 


* 

op. cit. 
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( 12 ) 


and 



LC C Q 

lap 



We 


(13) 


where the subscript s denotes the steady state assumption. 

Then the ratio of reaction rates predicted from the Inclusion of 
concentration fluctuations to those predicted neglecting these 
terms are 


Pc /Pt 

< d V dt >s 


= 1 + 


C 'C A 
a p 

C C ft 
a p 


(14) 


and an identical equation for the relative rates of depletion of 
the p* species. 


The limits on Equation (14) are readily determined from 
Equations (7) or (8) and elementary statistics. First_, we note 
from Equation (7) that for irreversible reactions, PC /Pt < 0 
and therefore u 


C 'CA 
a P 

C Co 
a P 


> 


Further, from elementary statistics we note that 

- 1 < 2_6 — < + l 




i — 


and therefore 


C ' 2 C ' 
a P 


2 \ b 


—2 —2 

c ci 

a P 


C 'C ' 

< < 


C ' 2 c ' 2 N 2 
a P 


C c B \ C 2 C 2 
a p ' a P 


(Cl) 


(C2) 


(C3) 


Substituting conditions (Cl) and (C3) into Equation (l4), we 
establish the limits 
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0 < 1 


ayat 

( a v 3t ) s “ 


1 + 





1 

2 


(C4) 


Conditions (c4) set the maximum errors in the prediction of 
reaction rates which the neglect of concentration fluctuations 
can produce. These limits are set by the individual variance - 
to mean-squared ratios of C a and Cg and are therefore func- 
tions of the marginal frequency distributions of C a and Cg . 


The limits established by conditions (C4) are shown graphic- 
ally in Figure 1 and we note immediately that the limiting errors 
in reaction predictions occasioned by neglect of concentration 

,2 


fluctuations are small for 



< 0.5 


but increase to 


highly significant values as this ratio exceeds 1.0. The poten- 
tial for order-of -magnitude errors in the prediction of the 
reaction rate exists whenever the product of the variance-to-mean- 
squared ratios greatly exceeds 1.0. 


The actual error depends, of course, on C'Ci/C a Cg . This 
actual error may be examined by forming the rat.io K 


r 


C a C P // C a 2 C J* 


(15) 


where r is the ordinary correlation coefficient and in this 
usage expresses the ratio of the actual error in reaction rate 
predictions to the maximum possible error for any given joint 
distribution of C a and Cg . Selected values of r are also 
graphed in Figure 1. In the^limit of r = 0 no error in reaction 
rate predictions is occasioned by the neglect of concentration 
fluctuations. This is, of cou rse , the situation when C and Cg 
are randomly distributed and C^Ci = 0 . However, it is clear P 
from Figure 1 that even modest values of r p roduce significant 

/ C f 2 C « 2 -\ 2 

errors in the reaction rate prediction when I — S E_ j s. i f 

C f3 / 


9 



to V 



Figure 1. The limits of errors of prediction 
of chemical reaction rates incurred 
by the neglect of concentration 
fluctuations. (See text for explana- 
tion of terms . ) 
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particularly when r < 0 . 

We shall return to. this analysis, and identify joint distri- 
butions of C a and Cg for which the fluctuations of concentra- 
tion must be included later. For the purpose of model develop- 
ment, we now turn attention to the importance of the third-order 
correlations in the chemical kinetics equations. 


The Role of the' Third-Order Correlation Terms 

Returning to Equations ( 9 ) to (11), it is evident that 
the primary role of the third-order correla tion s is to be found 
in their control of the rate of change of C 'Cl , both directly 

up p — p 

and through the rates of change of the variances, and Cg 

The effects of the third-order correlations on the reaction rates 
will therefore appear primarily as a time-integrated effect on 

2 

and any cumulative error in the estimates of and 

2 

C 'Cl will produce a cumulative error in C 'Cl . 
a p a p 

( 9 ) through (ll) 


(C5) 


their effect on C^C^ is negligible and we may close the model 

equations by setting C^Cl = C^Cd 2 = 0 . To illustrate that 
conditions (C5) are met under any^given circumstance, we must 
evaluate the joint distributions of C a and Cp from which these 
moments are derived since there are now no limiting conditions on 
their marginal distributions. In other words, we must turn 
attention to the distribution functions from which the means and 
moments have been derived if we are to determine the importance 
of third- or higher-order correlations in chemical reaction rates. 
Ideally, we would examine simultaneous experimental measurements 
of C a j_ and Cp^ to make this assessment; unfortunately, very 
few such data exist. However, so long as we assume that the basic 
chemical kinetic equations are correct (and our whole theory is 
based on this assumption) we can proceed by solving these basic 
equations for various initial distributions of C a; ! and Cg^ , 
determining in the process the time histories of all of the p 
relevant moments of these distributions. 


We may deduce immediately from Equations 
that if 


_ 2 

/"i n 1 c - 


and 


C ' Cl « CoC ' + C C 'Cl 

ap pa a a p 


— 2 — 

/-1 o 1 L 


C 'CA « C Cl" + C ft c 'Cl 
a p a p Pap. 
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The Moment Generating Model 


Under the assumption that only chemical reactions are opera- 
tive in changing the concentrations of a and (3 , the chemical 
kinetic equations can be written as total derivatives and 
integrated directly as a function of reaction time. 

dC . 
ai 

dt 

- - k l C al C PI 

( 16 ) 

dC P i 

dt 

- - k 2°ai C P i 

(17) 

and 



C . (t ) - C . (0) 
aa v ' ai v ' 

- Tt < c Pi (t) - C Pi (0)) 

( 18 ) 


while 


c Pi (t) 


c Pi<°> 


.( 0 ) 

k, ax v ' 


(19) 


1 - ^ 


k 0 C . (0) 

2 ai v y 


h exp 


( C Pi(°> - kf C cti^ ^ k f 


Equations (l8) and (19) specify the joint values of C a j_ and, 

Cpi at time t , given their initial values and the reaction; 
rate constants. They may be used to specify the frequency distri- 
bution of (Crd> CAj_ ) at any time t , given their initial 
frequency distribution, n-^(O) since, in the absence of mixing, 
n^ is conserved as C a i and Cg^ change value due to chemical 
reaction. Equations (lo) and (197 provide the information 
necessary to calculate all of the relevant moments of n i(C a i,Cg i ) 
and their rates of change. We may introduce any arbitrary initial 
distribution n i(C a q, Cp^) , subject only to the constraints 


and 


0 < C ai i 1 


0 < c P i < 1 


0 < ° a i + c pi < 1 




> 


J 


(C6) 
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In order to illustrate the general behavior of the distribu- 
tion of C a and Cp and the associated first-, second-, and 
third-order moments, we have chosen a simple distribution of 
points along the line C a = 1 - Cp and weighted each point 
equally (n^ = l/N) . The time history of n^C^i, Cpj_ ) is 
shown in Figure 2 and the moments of these distributions are 
plotted in Figure 3. Note particularly the distortion of the 
originally linear distribution of (C a , Cp) and the associated 
decrease (from zero ) of the t hird -order moments. The relative 

magnitudes of C^Cp and C pC^ + C" a C £ are al so plo tted in 

Figure 3 and, as can be seen the re , C^Cl and C 1C l 2 completely 
dominate the initial change of C 'Cl . However, this latter 
effect is too short-lived to be significant in the prediction of 
the time history of mean concentrations. This fact is shown in 
Figure 4, where the predictions of C a and Cp as a function of 
time, first, neglecting the fluctuations completely, and then 
neglecting only the third-order moments, are compared with the 
exact solution. The latter assumption produces an error of 
approximately 10 per cent at kt = 10 while the total neglect of 
the fluctuations produces an error of 300 per cent at that time. 

As a further example, and one which illustrates the impor- 
tance of the third-order correlations, we have constructed the 
distribution functions which simulate the case of intermittent 
sources. For physical perspective, imagine a free-way, oriented 
across the wind and on which the automobile traffic ranges from a 
steady, bumper-to-bumper stream to only an occasional vehicle. 

We assume that each vehicle emits approximately the same amount of 
pollutants per unit time, but that the ratio of the a and (3 
species emitted is slightly variable from one vehicle to another. 
Now we ask, "What is the average reaction rate for these exhaust 
materials immediately downwind from the roadway as a function of 
the intermittency of the traffic?" 

We simulate this situation by the frequency distribution for 
(C a , Cp) shown in Figure 3. The variability of C a and Cp due 
to variable exhaust emissions is portrayed as a circularly symmet- 
ric distribution and we take C a = Cp = 0 when there is no 
traffic upwind of our observation line. (Small background concen- 
trations have been assumed in another calculation but produce no 
significant effect.) We assume further that the pairs of nonzero 
values of (C a , Cp) occur with equal frequency and we measure the 
intermittency of the traffic as the fraction of time there is a 
vehicle upwind of the observation line, t^ . The moment -generat- 
ing model has been used to determine the chemical reaction rates 
as a function of kt for t^ = 1.0, 0.5, 0.33, 0.2, 0.1 (tq = 1.0 
corresponds to a steady, bumper-to-bumper stream of traffic). 
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Cot 


Figure 2. Example of the time history of the 
joint frequency distribution of 
(C a i, Cpj_) given the initial 
distribution shown for t Q and 
ki = k 2 = 1 . Each point was 
weighted equally (n^ =1) for 
calculation of the moments of 
these distributions (shown in 
Figure 3 ) . 
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Figure 4. Comparison of 
the predictions of C a 
and Cr , under assump- 
tions listed, with exact 
values from the moment 
generating model using 
the distributions shown 
in Figure 3 . 



Figure 5- Joint distri- 
butions of (C~, Cp) 
ahosen to simulate 
chemical reaction rates 
immediately downwind of 
a. roadway on which traffic 
is variably intermittent. 
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1 

2 


The values of 



r 


C 'CA 
a p 



^ 2 ) i 


and 


A~ = 


C ,2 C ' 
a P 


n n i 2 n 


C ft C ' 
p a 


+ C C^CT 
a a p 


at kt = 1 


and for each_value of t. are shown along with the observed 
ratios of c)C a /Pt to the reaction rates predicted assuming steady 
values of C a and Cg , (dC a /dt) s , and to the rates predicted 
neglecting the third-order correlations, (dCg/^t)j , in Table 1. 
These results, which are now firmly grounded in reality, show 
clearly that the potential errors in the prediction of reaction 
rates neglecting the concentration fluctuations can be realized. 
Although the effects of concentration fluctuations are negligible 
for a steady stream of traffic, tj_ = 1 , chemical kinetics based 
on this assumption underestimate the initial depletion rate of 
C a and Cp by a factor of nine when t^ = 0.1 


The effect of 
evident at kt = l 


neglecting the third-order moments is not . 

, however, since we started with the correct 


valu e of C 'CA and the integral effect of 

is nob yet_large at _kt = 1 . Their integral 


c i c 

pre 


2 

neglecting C ' CM 
effect on 


[ictions of C a and dCo/Pt can be seen, however, in 


and 
the 
the time 


histories of these quantities^ when the third-order correlations are 
neglected. Values of B-, = Co/C a j and B2 = (PC^/PtJ/^Co/Pt) j 


for each 
in Table 
sign. ] 


value of 

2 . [( 


J 1 . 
and 


kt 


) indicates 


_ a f 

C a l 


5, 10, ’16, an 
or (Sc c /St) I 


d 20 are listed 
has wrong 


An inspection of this table shows that the first effect of 
neglecting the third-order correlations, while retaining the first- 
and second -order moments, is to produce significant errors in the 
reac tion rate. This is, of course, the first integral effect on 
. The second effect is on the mean values predicted for the 
conc entrations , an error which depends upon the time integral of 
C^CM . We also note that not only can the errors be large (for 
example, the reaction rate is over-predicted by a factor of 100 
when tj_ = 0.1 and kt = 20) but_we get the totally erroneous __ 
results of positive values of dC a /dt and negative values of C a .' 

Prom these results we conclude that a generally useful 
chemical kinetic model must include the representation of the 
moments of the concentration fluctuations through the third-order. 
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Summary 


The various considerations and examples of the effects of 
concentration fluctuations on chemical reaction rates discussed 
above may be summarized as follows: 

1. The effects of concentration fluctuations can be signifi- 
cant , to the point of dominating chemical reaction rates. The 
situations under which they are significant are characterized by 
joint distributions of the reactant concentrations which are skewed 
toward large values of these concentrations , since this is the 
condition under which the variance-to-mean-squared ratios can be 
large with respect to one. 

2. The chemical reaction rate is significantly accelerated 
when the concentration fluctuations are positively correlated and 
is depressed when these fluctuations are negatively correlated. 
Combining these two requirements for significant effects of concen- 
tration fluctuations, we may sketch the general character of the 
joint distributions of C a and C g which require the inclusion 

of concentration fluctuations in chemical kin etic models. These 
are shown in Figure 6 . The first, when C^CJ > 0 is essentially 
the distribution of concentrations used for the intermitte nt 
traffic situation in the last section. The second, when C^Cl < 0, 
may be quickly identified with the situation when the chemical 
reaction is retarded or stopped by local depletion of one of the 
reactants. This latter case depends strongly upon the relative 
rates of chemical reaction and of local diffusive mixing. Further 
study of this situation depends upon an appropriate coupling of the 
chemical and diffusion equations. 

3. For strongly skewed distribution of C a and Cg , the 
case when concentration fluctuations become dominant in the deter- 
mination of the chemical reaction rate, the third-order correla- 
tions of these distributions must be included in the generalized 
chemical kinetic model. The neglect of these terms leads to the 
highly erroneous result that either the reactants are produced 
rather than depleted, or the mean concentrations go to negative 
values . 


CLOSURE OF THE CHEMICAL SUB-MODEL 


Having demonstrated the need for a more general chemical 
kinetic model for reactions in inhomogeneous mixtures, we may 
proceed along either of two fronts. The more general of these is 
to write down prediction equations for the third -order correlations' 
and seek closure by suitable assumptions regarding the fourth-order 
moments. If no such assumptions exist at the fourth order, we may 
proceed to higher-order moments until they are found. This line of 
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Figure 6. General 
types of joint 
frequency distribu- 
tions of C a and Cg 
for which concentra- 
tion fluctuations are 
significant in deter- 
mining chemical 
reaction rates. 
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inquiry has been pursued to the level of the fifth-order moments 
in the present study but has been set aside , primarily because the 
proliferation of simultaneous partial differential equations for 
the higher-order moments poses staggering computer requirements. 
This latter feature of higher-order closure becomes prohibitive 
when one considers the requirement for an equivalent level of 
closure of the diffusion equations with which this chemical sub- 
model is to be coupled. 


With these facts in mind, we have turned attention to the 
development of approximate closure schemes in which maximum 
possible information regarding the third-order correlations is 
sought from the first- and second-order moments. As is noted 
later, such a closure scheme ca nnot b e exac t. We begin by review- 
ing the definitions of C 'Cl , C ,2 C1 and C 'Cl 2 . 

a p ’ a p a p 

First, by definition, 


C 'C A = 
a p 


N 


n, (C 


ca 


" C a) 


(C 


Pi 


- V 


( 20 ) 


where nj_ is the frequency of occurrence of the joint values of 
(C ai , Cpj_) and N = 2 nj_ . Expanding Equation (20) and making 

use of the definitions of (L, and Ca , i.e., C. = i 2 n. C.. 
(j = a, P) , one obtains a P J N i 1 


C 'Cl 
a p 

C CV 
a p 


C C fi N i 
a p 


2 n, C ,C ft . - 1 
i ai pi 


( 21 ) 


We note in particular from Equation (21) that, since 


— - — 2 n. C ,Co, > 0 

C C ft N i 1 01 Pl 

a P 


( 07 ) 


C^CI/CqCq = - 1 only when the lower limit of conditions (07) is 
satisfied. This can occur only when any nonzero values of C * 
are coupled with zero values of Cg-? , and vice versa. Then when 
C 'Cl/C a Cp = - 1 , which is the condition for the termination of the 
chemical reaction, all joint moments of n i (C a j_, Cgj_) about the 
origin must be zero. 


Now, consider the definition of 




C ,2 C ' = 4 2 n. (C . 
a p N ^ i v ax 


- C c/( C Pi 


- c p) 


( 22 ) 
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C' 2 CA 
a p 

— 2— 

CTCn 
a p 


1 2 
— 2 n. Cr.Co. 

C 2 C r N i 1 ai Pl 
a p 


C i L C'CA 

1 + -2- + 2 -SL& 


C C R 
a p 


( 23 ) 


using Equation (20) and the definitions of C a and Cp . It is 
immediately evident from Equation (23) that our problem of approx - 

p 

imating C' Cd by first- and second-order moments reduces to find- 
ing a suitable approximation for the first term on the right-hand 
side of Equation (23) in terms of these moments. For convenience 
we denote this term by A. , i.e., 

JL CL 


A. = — ~ — 2 n. C 2 .C a . 

la C 2 C ft N i 1 a P . 
CL p 


(24) 


In seeking this approximation we may note the following conditions 

which A. must satisfy 
ia J 

1 . A . =0 when C~^Cj/C C ft = - 1 

ia a p a p 


when C ,2 /C 2 = 0 
a ' a 

c > 2 — o — 

-2- + 2 when C ,2 C ' = 0 

(T 2 C C A a P 

a a p 


4. A ± > 0 at all times. (C8) 

Condition 4 in (C8) operates primarily to constrain the joint 

values which /C a and CiC^/C^Cg may assume. For example, 

when conditions 3 and 4 in (Co)' are applied jointly, we derive 
directly the further condition that 
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(C9) 


2 2 

Finally, we must note that C ' Cl and C^Cl cannot be 
specified exactly from first- and Second-order moments alone. 
Consider, for example, the two distributions of C a j_ and Cgj_ 
shown in Figure 7. All of the first- and second -order moments are 
identical for these two distributions, but the third-order correla- 
tions are of opposite sign. Clearly, we are looking for useful 
approximations, not exact relationships, and any such approxima- 
tions must be delimited as to their regions of applicability. The 
criteria for usefulness of these approximations must be based on 
the accuracy with which the chemical model within which they are 
incorporated predicts the chemical reaction rates. However, in 
this development stage we shall be primarily interested in the 
accuracy of specification of the third-order correlations. 


Approximations of A^ a From First- and 
Second-Order Moments 

Utilizing conditions (C8) and (C9), it is not difficult_to 
formulate approxi mate exp ressi ons for A^ a and Aj_g from > 

Cg , C^Cl , and C^ 2 and C J 2 . For example, an early approxima- 
tion which has been tested is 


A 


ict 
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r 2 

C 'O' 


C f 

1, V 


a 

1 + _ 2 
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lo 

lo 


o 

L a P J 


a 


(25) 


Since C'Cl = 0 when C^ 2 = 0 and C^Cl/C^Cg ^ - 1 , it is 
evident from inspection of Equation f25) that p this approximation 
satisfies conditions 1, 2 , and 4 of (C8), but satisfies condition 


only if C 'Cl = 0 when C ' 2 C1 = 0 . According to condit ion (C9) 
the latterresult is admitted p but is not required, i.e., 

may be nonzero when C ,2 C1 = 0 . 

a p 


3 
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The development of an Improved approximation of Aj_ a 
improved in the sense of more appropriate realization of 
condition 3 in (C 8 ), proceeds as follows: Let 
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ia 




C l C ' 
2 a ft 

C C A 
a p 


2 2' 


f (C ,C A ,C 'C ' c ' c A ) 
cr or ?>’ a p' a p ' 


(26) 


where the function f satisfies the conditions 

a 


f = 0 when C <CA/C C ft = - 1 
a a p' a p 


f a = 1 when = 0 


(CIO) 


A simple function which satisfies condition (CIO) is 


f = 
a 


M + 1 


(27) 


where M is defined as the value of 0 'CA/C a Cg when C ' 2 CA = 0 . 
We may note immediately that M as defined must satisfy condition 
( C9 ) , i . e . , 


M > - \ 


As a further measure of M , however, we may also note that 

p P — 

since the ordinary correlation coefficient, C 'CA/(C 1 CAp 2 , must 
lie between ± 1 (condition (C 3 )) a p a p 
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(oil) 
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(C12) 


2 — 2 

and we expect M to be related to both C^ /C~ 
This expectation is reinforced by the fact that, 


2 

must have the same value for both C ' CA and 

2 — 2 U o ~2 

C' /C a is not necessarily equal to C^/C^ . 


and CA 2 /cI . 
as de fined, M 

C <^ 2 , although 
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Some values of M versus —%•— — , as determined from 

C Cq 

a p 

the moment generating models are shown in Figure 8. (Each of 
these solutions was derived from a log-normal or a composite of 
log-normal distributions of C ^ and Co. which were skewed 
toward large values of C a and Cp . ) p As can be seen from 
Figure 8, within the range tested, M is a relatively well- 

r — - — - n , 


behaved function of 


2 2 
C i* CJT- 
a p 

—2 —2 

C^ C R 
CL p 


. For our present stage of 


approximation, however, we have chosen a dichotomous relationship 
for M , namely, r 


M = 0 when 


M = 1 when 
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?j 2 tt 2 

°a °P 


2 2 
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C C R 
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(C13) 


Using Equations (23), (26), and (27), our approximate predic- 

2 ? 

tor equations for C 1 Cl and C 'Cl are 

a p a p 


c ,2 _ C a°P 

a P p “ M + 1 
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C 'Cl C 'Cl 
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and we specify M by conditions (C13). 


C Co 
a P 


C 'Cl 
a P 

C Co 
a P 


( 28 ) 


(29) 


As a first, but severe test of this approx imation , comparison 

of the predicted and observed values of C^C^ for the most 

2 p i 
C ' C XT 2 

extreme value of shown in Figure 8 (Run L24) is shown 

c c* 

a P 
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Figure 8. Values of M vs — 

C 2 C 2 
a p 

as determined by solution of the 
basic chemical kinetics equations 
for various initial log-normal 
frequency distributions of C„ 
and C Q . 
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in Figure 9> and the predictions of PCo/Pt from this model are 
compared with those predicted by 1) neglecting the fluctuations 
completely and 2) neglecting only the third-order correlations, in 
Figure 10. 


The Approximate Chemical Sub-Model 
For Inhomogeneous Mixtures 

Employing Equations (28) and (29) and Conditions (C13), 
the chemical sub-model for two-body reactions in inhomogeneous 
mixtures is 
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(31) 
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(33) 

( 34 ) 
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Figure 10. Comparison of predicted vs actual chemical 

reaction rates for Run L24 under the follow- 
ing assumptions: 

1. Neglect of concentration fluctuations 
(uniformly mixed) 

2. Neglect of third-order correlations 
(Mark I) 

3. Inclusion of approximate estimates of 
third-order correlations (Mark VII ) 
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( 36 ) 


(C14) 


and the sub-model Is closed at the level of the first- and second 
order moments. 

The only rigorous test of this approximate model available to 
us now is a comparison with the exact solutions of the chemical 
kinetic equations, as a function of reaction time kt , and for 
various initial distributions of (C a j_, Cpj_) . The much more 
realistic case of coupled chemical depletion and dilution by turbu- 
lent mixing must await the coupling of the chemical and diffusion 
sub-models. However, if the local diffusive mixing is very, very 
slow compared with chemical depletion, this sub-model must "track" 
the chemical depletion correctly. The following tests of the 
chemical sub-model are therefore restricted to this circumstance 
so far as any degree of reality .is concerned. 

Since our primary concern in the coupled chemistry-diffusion 
models will be accuracy in the prediction of the local mean 
reaction rates, we are particularly concerned with this facet of 
the chemistry sub-model. The comparison of Interest is between 
the local depletion rates of the reacting chemical species as 
measured by bc a /bt . From a variety of initial distributions for 
n i( c ai' c Si) t we have chosen four which exhibit varying degrees 
or the effects of inhomogeneous mixing. Their initial distribu- 
tions are shown in Figure 11 and the comparisons of reaction rates 
as a function of reaction time are shown in Figures 12 to 15, 
(including the initial reaction rates predicted when the concentra- 
tion fluctuations are neglected). 

These comparisons, although by no means exhaustive, show that 
the approximate chemical sub-model developed here captures a very 
large fraction of the effects of inhomogeneous mixing on chemical 
reaction rates. Over a very wide range of reaction rates, this 
model predicts the exact rate to within a factor of two, while the 
neglect of the fluctuation terms in the chemical kinetic equations 
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Figure 12. Comparison of reaction rates 
predicted by the chemical 
kinetic model (Mark VII ) with 
exact solution for initial 
distribution of (C n , Ca) shown 
in Figure 11. P 
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Figure 13. Comparison of reaction rates 

predicted by the chemical kinetic 
model (Mark VII ) with exact solu- 
tion for initial distribution of 
(C , Co) shown in Figure 11. 


34 




Figure 14. Comparison of reaction rates 

predicted by the chemical kinetic 
model (Mark VII ) with exact solu- 
tion for initial distribution 
shown in Figure 11. 


35 



aC a /at 



Figure 15 . Comparison of reaction rates predicted 

by the chemical kinetic model (Mark VII) 
with exact solution for initial distri- 
bution shown in Figure 11. 



incurs of up to a factor of 100. On this basis it seems safe to 
proceed to the coupling of this chemical sub-model and the invari- 
ant (second-order closure) diffusion sub-model. 


THE CONSTRUCTION OF A TWO-DIMENSIONAL 
COUPLED DIFFUSION/CHEMISTRY MODEL 
FOR A BINARY REACTION SYSTEM 


The development of the closure scheme for the chemical sub- 
model described in the previous section, coupled with the models 
for prediction of turbulence structure, fluxes, and turbulent 
diffusion of matter provides the necessary modules for a coupled 
diffusion/chemistry model, the objective of this program. In 
particular, these developments make possible coupled models which 
permit examination in detail of the processes of generation and 
decay of the fluctuating components, as well as the mean values, 
of turbulent diffusion and inhomogeneous reactions of a binary or 
two-body reactive system. 

As a starting point for these coupled models we have chosen 
the relatively simple, but realistic, situation of two reactive 
but otherwise passive pollutants emanating from either common or 
separate cross-wind line sources. This .choice reduces the diffu- 
sion calculations to two dimensions and permits the decoupling of 
the diffusion/chemistry model from the turbulence model, since 
there is no feedback into the turbulence field due to either pollu- 
tant density or exo- or endothermic reactions. The decoupled 
turbulence model is used tO|, generate the field of turbulent motions 
and fluxes, which are then 'jused, along with 1 ; source specifications 
and reaction rate constants, as input to the coupled diffusion/ 
chemistry model. As can be seen from the full derivation of the 
two-dimensional model presented in the next section, this system, 
involving as it does nine simultaneous partial differential equa- 
tions, is already rather complicated. However, one of the primary 
reasons for starting at this basic level of complexity is to permit 
close examination of the interactions between turbulent diffusion 
and chemical reactions in the simplest realistic mode in which they 
could occur. When these interactions are understood, extension of 
the models to three-dimensional configurations, nonpassive pollu- 
tants, and three -body reactions can be undertaken with a much 
better appreciation of the complex nonlinear system within which 
they will operate. 

Due to the limited time and funds for this project, only a 
few test calculations of the combined effects of turbulent diffu- 
sion and chemistry have been possible. Some basic calculations, 
which begin to define the effects of turbulence vis a vis chemical 

reaction rates, and two sets of calculations for the NO - 0 

x 3 
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patterns to be expected in an SST exhaust plume are presented and 
discussed later. Despite their limited number, these examples 
already point up sharply the effects of inhomogeneous mixing 
(produced by the turbulence field) and the effects of diffusion- 
limited conditions on fast reactions characteristic of photo- 
chemical chains. 


Derivation of the Modeled Equations for the Mixing of 
Two Chemically Reacting Materials Emanating 
from Cross-Wind Line Sources 


For an atmospheric shear layer in which the Schmidt number is 
equal to one and the adiabatic density is constant, we may follow 
Donaldson* and write the equation governing the diffusion and 
chemistry of a reacting species, a , as 


SC 

a 

SV 


SC 

u j 3x7 + v o 

J 


a 2 c 

S 


a 


k C C ft 
a a p 


(37) 


where k a is the reaction rate of the a species with a second 
species S . A similar' equation may also be written for Cp , but 
this will not be done until the modeling is completed. ! 

We may express our variables in terms of their mean and fluc- 
tuating parts as 
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+ u ! 

3 


(38) 


Substituting these expressions into Equation (37) and averaging, 
we obtain the mean local rate of change of concentration in terms 
of convection, molecular and turbulent diffusion and chemical 
reaction 
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In deriving Equation (39) we have used the continuity equation 
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(4o) 


From Equati on (3 9) we see that we now require expressions for 


u !C^ and 


C'Cl . We first subtract Equation (39) from Equation 

~ ~ ^ -t- ' ^ ^ 4-1 -O “I . , ~ 4 ,->4-4 /If 


( 37 ) to obtain an equation for the fluctuation , 
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(41) 


Following Donaldson , we may write the expression for u^ as 
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(42) 


Multiplying Equation (42) by C^ and (4 l) by u| , adding, and 
averaging, we obtain an expression for u.'C ' , 

1 Qj 
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Equation (43) introduces various third-order correlations, pressure 
correlations and dissipation terms that have bee n mode led previous- 
ly. We must, however, obtain an expression for C^T 1 . We may 
write the equation governing the temperature fluctuation T 1 as 
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We then multiply this equation by , Equation (4l) by T* , add 

them together and a verage to obtain the equation governing the rate 
of change of CH? 1 
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The triple correlation 
already been modeled, 
below . 


u-jC^T 1 and the dissipation term have 

The correlation T'C'Ci will be discussed 
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Returning to Equati on ( 39)* we see that we must determine the 

governing equation for . If we multiply Equation (4l) by 

Cg and add to it the equation obtained by multiplying by the 

fluctuation equation for Cp ( a replaced by P and p by a 

in (4l) ) then the average of that expression gives the equation 

for C 'C ' 
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We see from Equation (46) that we must model the third-order 

2 2 

correlations C^C^ and , but also that we require expres- 

expressions for and . These two equations a re i denti- 

cal except for the transposition of a and p . The C' 2 
equation is found by multiplying (4l) by and averaging to 

obtain ' 
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Now, since we are dealing with an atmospheric shear layer, and 
our initial attention will be directed to cross-wind line sources 
of a and p , we expect that the_ only derivative of importance is 
the one normal to the mean flow u in the x cartesian direction. 


4l 



Thus, only x j = x, = z will he important in the equations. Also, 
we may set t = x/u without loss of generality. The modeling 

for the third-order correlations, pressure correlations, and 
dissipation terms is prescribed by Donaldson as 
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(48) 


and X, A]_, A2 and A^ are 

a 


where q = (u'u 1 + v'v' + w'w') 2 

length scales. Also, A c = (Ao Ao ) 2 where a and (3 

2 B 2 cx 2 p 

correspond to the species in question. When Equations (48) are 
substituted into Equations ( 39 ), (43), (45), (46), and (47), we 
obtain our final equation set governing the simultaneous diffusion 
and chemical reaction of a and p . 


u 


PC 


b 2 C 

a a 

= V °^ 


bC ' 


a 


w 




- k. 


n (C C ft + C 'Cl) 
1 v a p a p 7 


(49) 


- Sc p 
U = V 


3% 


P 


bz l 


bCAw 1 
P 

3z 


- k, 


2 ( C a C 6 + 


(50) 


42 




* *S5T 


_ . - T 

_ VJ ,VI oZ 


a + '■'a, 

o 


+ A, ) * ^ 


b ((,2^2 

+ V 


a 


} a 


- 2v 


q. ■ tr T ^ 7 + ° 3z 

Ct 

H 


c a*l 

a 


S^flX 15 ^ 


X 


1 


C a C r 


_ fi, <i "5z 

CaP^ 2 B 


t + ^gra 


„ ^vT’ Sz o 


u 


ip 


*c§* 


f fv ^ ^ "§£ 

* + 3 %^r ^£ 

_ , 2v ^ 


O \ 


CL- Tfe? 7 + 0 3z c 


3 


1\ 


1 


'P 


$ 


^ak\ ^ 2) 


. K,\*<ft? 


T + 


T - ho ^ oZ 
^3 


&T 


i^%p _ - “ST 


~\fn 1 /Vn ^ **5^ 

^ c a „ If^ 7 ~Sz + 


'a 


a 


u -Sx 


a 

+ v o 3z 2 


3 2 ^f - 2v 


O X 2 


a 


(. 53 ) 


. K, 


{*& 


- . rfi^Jpp, 

jr -nT^ + a P 

TTJT + a. 


^3 



dcTF 1 


ac 


p 


a? a 


- ' ^ <£* - W & + & 1% q 3* 


acjr 


a C'T' C'T' 

+ v„ —I# 2v -V 


a 2 ‘ 


x 


0 


- ko < c c7T t + + T 'C 'C ' 

^lap pa ap 


(54) 


acHcj 

“33^ 


t acp 


ac 


C’w 1 -3-C- - C'W‘ *t— — + 4- 

a oz P oz dz 


ac 'C ' 

JD 


+ V 


&2c 7p 
° ip 


- 2v 


C 'C ' 
a p 


o 


X 


J 


- _ .2 - 


k, C C ' + C ft C 'C ' + C 'C ' 

1 l a P P a P ap 


- „ .2 


kp 1 C C 'C ' + C R C • + C ' C ' 

^ l a a P pa ap 


(55) 


ac 


— a 

u 33^ 


ac 


ac ,£ 

a 


23Tw T ^ ^ ) A 2 q 3z 


a 


3 2 C ,2 C' s 

+ v -t£- - 2v -2, 

° 3z 2 ° 

a 


a 


~ , 2 


- 2k-, < C C 1CX + C A C ' + C ' C ' 

llaap pa a P 


(56) 


44 



- ^C ' dC R 

u 33T-= - 2C ^' - 


p w sr* + 3F 



(3 C 


+ V. 




° dz' 


- 2v 


5L 

° X 2 

P 


- 2k, 


_ 2 — 


C Cl 
a p 


c ft c -ci 

P a P 


+ C 


•C 

a 


i 


X 

j 


(57) 


We have written k a s k-^ and kg = k2 . Since we have decoupled 
the diffusion/chemistry model from the background turbulence model, 
a solution of Equations (49) - (57)_requires knowledge of the 

initial d istri butions of C a and Cg and of the flow parameters 

u, T, T 0 , w'w', q, and w'T 1 . The macro-scale lengths A and 
micro lengths X must also be known as functions of the background 
turbulence or the plume characteristics. 


Finally, the coupled diffusion/chemistry model is closed by 
modeling the third-order chemistry correlations as described previ- 
ously. 
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A gen eraliza tion of (58) and (59) gives the appropriate modeling 

for T 1 C 'Cl 
a p 
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With the properly modeled equations, we can now proceed to a 
discussion of some of the results of computer solutions of these 
equations . 


SOME CALCULATIONS OF THE INTERACTIONS 
OF TURBULENT DIFFUSION AND CHEMISTRY 


As must he evident from the derivation and recapitulation of 
the closed diffusion/chemistry equations, it is virtually impossi- 
ble to trace the effects of variations of any one variable through 
this simulation system. This being the case, the validation of the 
model must involve multiple iterations with a systematic, step-by- 
step variation of each of the input variables, and comparison of 
these model predictions with observed values of concentration 
patterns, chemical depletion rates, turbulent flux divergences, 
and the like. Neither time nor resources have permitted this 
validation of the model, of course, and the calculations presented 
here must be regarded as suggestive rather than authoritative as 
to how real chemically reactive and turbulent systems may operate. 
The verification of these predictions must be deferred, but the 
results which are presented here argue strongly that at least in 
some circumstances the interactions of turbulent diffusion and 
chemical reactions are highly significant and, if verified from 
observations, models of this type may also improve predictions of 
air quality in the lower atmosphere significantly. 


An Illustrative Calculation 

In view of novelty of simultaneous consideration of the 
turbulent diffusion of reactive chemical species and their 
reactivity in inhomogeneous mixtures, it appears highly desirable 
to examine in detail the individual processes by which turbulent 
diffusion and chemical reactions produce observed patterns of 
reactant concentrations and reactant depletion in a simple but 
realistic system. To this end we have chosen to calculate the 
combined processes of diffusion and chemistry for the case of a 
plane jet of reactant a released continuously and isokinetically 
into a uniform environment of reactant (3 . The environment of 
|3 is characterized by a uniform transport speed u = 10- m/sec 



and an isotropic and homogeneous field of turbulence characterized 

p 2 2 

by the vertical intensity of turbulence w 1 = 1 m /sec . The 

plane jet of the a species is oriented across the mean field of 
flow and the initial vertical distribution of the concentration of 
a, C , is taken as gaussian with a central value of one and stand- 
ard deviation = 0.4 m . In view of the requirements that the 
mass fractions of a and (3 equal one, this "jet " of the a 
species displaces the ambient (3 species at the source in such a 
way that the initial distribution of the_concentration of the (3 
species is the complementary gaussian, Cp 0 = 1- C nn . This geometry 
of the initial distributions of C a and Cp is shown in Figure 16. 
Note that no initial fluctuations of C a and Cp are introduced 
at the source. Finally, we take k^ = ^ ~ 1*0 • 

In keeping with the constraints imposed in the construction 
of the model, we assume the reaction of a with (3 proceeds 
isothermally . For our present purposes we shall also assume that 
this reaction is irreversible, even though this assumption is not 
mandatory*. With these input conditions the model calculates the 
redistribution and the chemical depletion of a and p as a 
function of travel distance or time after emission. 

As a first partial view of the coupled effects of diffusion 
and chemistry in this flow reactor, we may compare the predicted 
axial concentrations of the a species as a function of distance 
from the source and under the following conditions: 

1. a does not react with |3 (diffusion only) 

2. Diffusion and chemistry occur, but the chemical 
reaction rates are calculated on the basis of 
the local mean values of concentration only. 

(We hav e ter med this "homogeneous chemistry " 
since is neglected. ) 

3. Diffusion and inhomogeneous chemistry are 
operative . 

This comparison is sh own i n Figure 17 and it is immediately evident 
that the neglect of leads to a significant over prediction 

of the rate of decrease of the axial concentration of the a 
species. For example, the ratio of the predicted concentrations 
at x = 40 m is two, and, as is clear from Figure 17, this ratio 


The effects of reversibility of reactions and catalytic cycles 
may be accommodated to a certain extent by appropriate choices 
of the reaction rate constants, k]_ and k2 . Similarly, 
three body reactions of the type - k^C a CpC]yj may also be 
simulated, if dC M /dt ~ 0 , by taking k = kj_C M , where M is 
the third body. 
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Figure l6. Source configuration for a plane jet 

of pure a species injected isokinetic- 
ally into an environment of pure (3 
species . 




Figure 17 . Comparison of the axial concentrations 
of the a species as a function of 
travel distance from the plane jet as 
estimated for 1) diffusion only, 

2) inhomogeneous chemistry plus 
diffusion, and 3) homogeneous chemistry 
and diffusion. (See text for details.) 
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is i ncrea sing with x . This result reflects primarily the effect 
of C{C1 on the chemical reaction rate. However, such a simple 
portrayal of the results of coupled chemistry and diffusion does 
not portray the balance of the diffusive and chemical processes at 
work. In order to gain this insight we must examine in detail the 
balance of turbulent diffusion and chemical reactions going on 
across the plume. 


In order to examine this balance we have neglected the 
molecular diffusion terms_and plotted each of the rates which 
determine dC^/Pt and dCp/3t as a function of distance from 
the plume centerline at x = 37 m or t = 3.7 sec. The profiles 
of mean concentrations of a and p are shown in Figure 18 and 
the diffusion and chemical reaction rates predicted by the model 
are shown in Figure 19. In order to discuss and interpret these 
results, we recall the balance equations for 3c a /3t and dCp/3t . 
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as small in 
Each of the retained 


terms is plotted as a function of height above {or below) the 
plume center line in Figure 19. (Recall that 3Cp/3t and 3Cp/3t 
are the local rates of change of the mean concentration of a an 


P due to both turbulent diffusion and chemical reaction; 

3 


and 

1— C~w ' 
oz a 


and 37 C p w 

the 


are the local divergencies of the turbulent flux of 

a and P species; - k(C’ a C’g) is the average local chemical 
reac tion rate due to the local mean concentrations of a and p 
- k C 'Cl is the average local chemical reaction rate due to 
correlated fluctuations of the concentrations of a and P . ) 


A patient inspection of Figure 19 reveals the following facts 
regarding the diffusion and chemistry processes through the plume: 

1. The diffusion of the a species is removing a from the 
plume core from the center line to z = 1.25 m and is causing an 
accumulation of a from 1.25 to 5 m (Curve 4). 

2. The diffusion of the p species into the plume is oper- 
ating to remove p from the height zone 2 to 8 m and accumulate 
P in the height zone 0 to 2 m (Curve 6). 
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Height above plume 







Rotes for terms indicated (ppp/sec) (x I0 2 ) 


Figure 19. Calculated values of local rates of change 
of the concentrations of a and |3 due 
to diffusive flux divergence and to 
inhomogeneous chemical reactions. 
(Concentrations are expressed in parts 
per part = ppp.) 
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3. The chemical reaction is depleting both the a and f3 
species between the plume center line and z = 6 m (Curve l), the 
upper height being the limit of a penetration into the p 
environment at this time. 

4. The chemical reaction diminishes the rate of increase of 
the concentration of P below z = 1.5 m and accelerates the 
decrease of concentration of p from 1.5 to 5 m (Curve 5). 

5. Below 2.75 tn the chemical reaction accelerates the deple- 
tion of a in the plume core, but above 2.75 m the reaction rate 
very nearly balances the diffusive transfer rate for a , i.e., 
above z = 2.75 m the chemical reaction is diffusion limited 
(Curve 3 ) . 

6. Below z = 5 m the diffusive mixing of P with a becomes 
increasingly inhomogeneous so that at the plume center line the 
chemical reaction rate is proceeding at only 60 per cent of the 
rate computed on the basis of mean values of P and a concentra- 
tion at that height. (Comparison of Curves 1 and 2) 

From these detailed comparisons we can immediately deduce that 
the a plume is not only being rapidly depleted by reaction with 
P but that it is also growing in vertical width only slowly 
because of the balance between diffusion and reaction rates in the 
outer limits of the plume. On the other hand, the P deficit in 
the initial plume is being filled in with only minor interferences 
from the local chemical reaction with a . Further, the depletion 
of the a species in the core of the plume is significantly slower 
than would have been expected from mean-value chemical kinetics. 

In this region the chemistry is "mixedness " limited while above 
this region it is clearly "diffusion" limited. 

This simple example is intended only to Illustrate the balance 
of diffusive and chemical processes (and, incidentally, points to a 
laboratory experiment which may verify these predictions). However, 
the example does illustrate the model 's power to simulate and 
portray complex chemistry/diffusion processes. 


The Sensitivity of Chemical Reactions 
to Turbulent Diffusion Rates 

The example discussed above is, of course, only one particular 
case and cannot provide any real insight into the sensitivity of 
the combined diffusion and chemical depletion of reactive species 
to variations in the input variables. Extensive sensitivity 
analyses have not been possible, but we have done a partial analy- 
sis of the effect of turbulence intensity on the chemical depletion 
rate of the a species for the flow reactor described above. 

Returning to Equation (62) we note that since u is constant 
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.is the horizontal flux of the a species at distance 
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and since the first term on the right of Equation (65) is zero, 
the rate of change of the flux of the a species is determined 
by the total reaction rate over the height of the plume at any 
distance x . In order to compare the model 's predictio n of this 
rate against a_ limiting condition, we may note that as C 'CA goes 
to zero and Co tends to be uniformly distributed through the a 
plume (a condition which can only be approached asymptotically), 
the basic_chemical reaction goes over to a first-order reaction 
and for u = constant 


1 



F dx - 
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( 66 ) 


where Cn Q is the environmental concentration of the p species. 
We can rewrite Equation (65) as 
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and compare Equations (66) and (67) from the model calculations 

2 

using various values of w 1 


This comparison is shown for the flow reactor problem and at 
x = 40 m in Figure 20. As can be seen from Figure 20, the relative 
chemical depletion rate of the a species is quite sensitive to 

2 

the intensity of turbulence for small values of w ' but becomes 

2 

quite insensitive to this input parameter when w 1 becomes 
larger. We may also note that even with the vigorous turbulence 
2 2 2 

of w 1 = 3 m /sec , the chemical depletion rate has only achieved 
about 60 per cent of the limiting, first-order reaction rate at 
x = 40 m or kt = 4 sec . 
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w /2 (m 2 /sec 2 ) 


Figure 20. The partial dependence of the chemical 
depletion of the a species on the 
intensity of turbulence. 



An Application of the Model to the Case 
of the Stratospheric SST Jet Wake 

With the successful development of the first version of a 
coupled invariant diffusion and chemistry model, we are in a 
position to run trial calculations of the combined effects of 
turbulent diffusion and chemical reactions as they may effect the 
distribution and chemical alteration of jet engine exhaust products 
in the stratosphere. In view of the fact that there has been no 
opportunity to validate the model against reliable experimental 
data, these initial calculations must be considered suggestive 
rather than authoritative descriptions of how stratospheric jet 
wakes behave. This caveat is doubly reinforced by the limited 
number of cases which have been calculated to date. (An extensive 
sensitivity analysis and cataloging of possible ranges of strato- 
spheric conditions is proposed for the next phase of this work. ) 
Definitive results must await these further calculations and, hope- 
fully, experimental verification. Nevertheless, even these early 
and partial calculations have already brought to light the very 
real possibility of interactive effects between diffusion and 
chemistry which significantly alter prior estimates of the dura- 
tion and chemical makeup of the SST jet wake. In particular, the 
following calculations begin to portray an intact exhaust plume 
which may persist for substantial periods of time and which effect- 
ively shields itself from highly reactive environmental constituents 
such as 0^ . 


Estimates of the Turbulence Fields 
in the Lower Stratosphere 

Since an important input to the coupled diffusion and chemis- 
try model is the field of turbulence and turbulent fluxes of heat 
and momentum, a series of calculations of the equilibrium values 

2 2 2 

of th e turb ulen t energies, u' , v 1 , and w 1 , and the fluxes 
w 'T 1 and u 'w 1 , using Donaldson's invariant model* were completed. 
A mean temperature lapse rate of + l°C/km was assumed and the 
mean wind shear was varied from 5 x 10 - 3 sec - -*- to 5 x 10 2 sec - -*-. 
(For wind shears less than ~ 5 x 10 - 3 sec -1 the critical Richardson 
number is exceeded and no equilibrium turbulence can be maintained.) 

The final input to the turbulence model is the integral scale 
length. For the initial calculations this length was taken as 15 m, 
primarily on the basis of verifiable results in the lower atmosphere 

2 2 2 

The values of u' , v 1 , and w' which associate with this scale 
length and with various values of wind shear are shown in Figure 21. 
As can be seen there, the equilibrium turbulence is strongly 


op. cit. 
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T urbulence intensity , 



Figure 21. Dependence of turbulent kinetic 
energies on wind shear at 20 km 
altitude for a temperature lapse 
rate of + l°C/km and a macro 
scale length of 15 m. 


dependent upon the mean wind shear, but, for reasonable values of 
shear in excess of that required to maintain turbulence,, the 
turbulence intensity is of the order of 10 3 to 10“2 m^/sec^. 

In order to test the sensitivity of these results to the 
integral scale length, the turbulence calculations were rerun with 
a temperature lapse rate of + l°C/km , wind shear equal to 10 2 
sec _ l (a rather large value), and the scale length equal to 100 m. 

p p p 

The equilibrium values of u 1 , v ' and w 1 for these conditions 
are shown in the following table. 


Scale Length 
m 

,2 
u ' 

2 / 2 
m /sec 

,2 
v ' 

m 2 /sec 2 

w 2 

m 2 /sec 2 

15 

-2 

1 x 10 

5 x 10‘ 3 

5 x 10~ 3 

100 

4 x 10" 1 

2 x 10 _1 

2 x 10 _1 


From these turbulence calculations we have chosen 

w 1 2 = 0.02 m 2 /sec 2 , w 1 T 10 ^ °C/sec, and u 'w ' = 10 2 m 2 /sec 2 
as representative of the turbulence field in the 18-20 km height 
range of the stratosphere. In a separate calculation, this choice 
was found to correspond to an effective eddy diffusivity, k zz of 
5 x 102 cm2/sec, a number which agrees well with Lilly's calcula- 
tions of persistent diffusivities as derived from Air Force HICAT 
data . However, for many periods of time these values of turbu- 
lence and diffusivity may be excessive, and the consequences of 
persistent, near-laminar flow regimes must be investigated, 
particularly since the effects of mixedness and limited diffusive 
transfer on chemical reactions between plume and environmental 
reactants are enhanced by reduced levels of turbulence. 

The further study of the effects of persistent turbulence 
fields and diffusivities is the subject of Contract NAS1-118'73 j 
with the choice of a representative turbulence field mentioned 
above, we may now turn attention to the simulation of an SST 
exhaust plume during the first hour or so after the aircraft 
induced vortical motions have subsided. The behavior of the plume 
during this period is particularly pertinent to the wake sampling 
of exhaust products which has been planned by NASA and a task to 
which the present contract work has been supportive. 


D. K. Lilly, Draft of D0T/CIAP Monograph 3. (Available from A.R.A.P. To be 
published in SIAM J. ) 
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The Simulation of Combined Chemistry and Diffusion 
in the SST Wake at Cruise Altitude 


As a first application of the two-dimensional diffusion/ 
chemistry model to the SST wake problem, we have chosen a simpli- 
fied version of the N0 X catalytic cycle for the destruction of Oo 
and, simultaneously, a calculation of the vertical diffusion of an 
inert exhaust constituent such as CO . We assume the aircraft 
flies along the direction of the mean wind at 20 km altitude in 
the sunlit sky. Then, assuming that the production rate of 0^ 
in the exhaust plume during the early post -vortex stage is 
negligibly small compared with 0o destruction by reaction with NO, 
and that the photolysis of NO2 ana the abundance of 0 in the plume 
are sufficient to restore all N0 2 to NO + Og, we may simplify the 
N0 X - O3 catalytic cycle to 

NO + 0 3 N0 2 + 0 2 (68) 


N0 2 + 0 -*■ NO + 0 2 


and 


S[0 3 ] 

3t 


chem ^1 


[0 3 ] [NO] + [0 3 J. [NO]. 


( 69 ) 

(70) 
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chem 
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The reaction rate constant, k^ , is taken as 5 x 10^ (1/ppp-sec). 


With this choice of the chemical reactions, NO, of course, 
serves as a surrogate tracer for any nonreactive exhaust product 
since the total amount of NO in the plume is conserved. We can 
therefore follow the history of the distribution of NO concentra- 
tion in order to determine the effects of turbulent diffusion alone. 


In order to establish the initial distributions of O3 and NO 
in and around the plume, we assume the NO is distributed initially 
in a gaussian profile about the axis of the plume, with a standard 
deviation of 25 m and a central value of 2 x 10"°ppp. We further 
assume that the combination of entrainment and photodissociation 
of NO2 during the vortex phase results in the 0o being initially 
distributed in a complementary gaussian profile with an axial _r 
value of zero and an ambient value, outside the plume, of 2 x 10“° 
ppp. Both of these initial concentration distributions are taken 
as steady mean values, i.e., all concentration fluctuations are 
initially zero. Finally, we take the density of the atmosphere as 
10 - 1 kg/m3 and the kinematic viscosity of air as eight times its 
sea level value. 
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With these inputs, the coupled model has been run to simulate 
the diffusion of NO and the diffusion and chemical depletion of 0.-, 
for a real-time period of 3000 sec. We begin the analyses of these 
results by examining the details of the calculation at t = 150 sec 
as they are portrayed in Figure 22, by plots of the concentration 
profiles of 0o and NO and the profile of local depletion rates of 
0o, and in Figure 23, where the profiles of the turbulent diffusive 
fluxes of NO and 0o and of the second-order correlations of 
concentration fluctuations are shown. 

Beginning with the mean concentration profiles in Figure 22, 
we see immediately that during the first 156 seconds the chemical 
reaction of NO with 0o in the plume has seriously depleted the O3 
in that region; the exhaust plume is essentially O^-free to a 
height of about 50 m above (or below) the plume center line; this 
was true only at the plume axis initially. Clearly, during this 
initial period the chemical reaction has been limited by the rate 
of diffusion of 0o into the plume and, as is shown on the right- 
hand side of Figure 22, this supply rate of O3 to the NO is such 
that all of the reaction is occurring in the outer fringes of the 
NO plume. In this case the limitations of diffusive transfer have 
operated to create a significant 0o hole, and the bulk of the 
exhaust products are contained in an O^-free region. From the NO 
profile we can see that the exhaust products have been diluted by 
about 10 per cent by the entrainment of environmental air. 

A comparison of the chemical reaction rates as estimated, 
first by only the mean values of the concentrations of NO and 0o, 
and then by inclusion of the "mixedness " term, > shows that 

in this early stage the reaction of 0o with NO is also mixedness 
limited. The ratio of the total reaction rate (integrated over 
the plume height) to that estimated from mean-value kinetics alone 
is 0.39; i.e., at this time the reaction is proceeding at about 
40 per cent of the rate which should obtain if the O3 and NO were 
well mixed locally. As is evident from Figures 22 and 23, this 
effect of "folding " of the two reactants with only partial local 
mixing is associated primarily with the plume edge, where NO is 
diffusing into the ambient On field and with the maximum flux of 
O3 into the NO plume. We can also note at this point that the 
diffusion limitation and the mixedness limitation are operating 
against one another; if the mixedness limitation were not present 
the Oo-hole would be even more prominent due to the more rapid 
depletion of O3 in the outer fringes of the plume. 

With this "snapshot " of the state of affairs at t = 150 sec 
in hand, we now proceed to t = 3000 sec. The profiles of NO and 
0^ concentrations and the chemical reaction rates at this time are 
displayed in Figure 24. We can note immediately that the diffusion 
limitation on the chemical reaction of NO with O3 is still opera- 
tive; the O^-hole is now evident out to 180 m from the plume center 
line and the chemical reaction is still confined to the outer edges 
of the NO plume, even though it has now been diluted to about 
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Height above plume , z (meters) 



Figure 22. Vertical profiles of the mean concentrations 
of O 3 and NO and of their local reaction 
rate at 150 sec after vortex decay of 
simulated SST exhaust plume. 


61 



Figure 23. Vertical fluxes of Oo and NO an d of the 
chemical mixedness term, C'CA , for 
the situation described in Figure 22. 



Height above plume (£_ , i (meters) 



Figure 24. Vertical profiles of 0 3 and NO 

concentration and local chemical 
reaction rate for simulated SST 
plume at t = 3 OOO sec . 
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4/10ths of its original concentration on the plume center line. 


The Sensitivity of the Oo-Hole to Variations in the 
Amount of NO in the Plume and to the Intensity of 

Turbulence 


Recognizing that the above calculation represents only one 
possible combination of conditions, we have completed a series of 
calculations aimed specifically at determining the combinations of 
diffusive transport of O3 and the amount of NO in the plume which 
permit the environmental Oo to penetrate to the plume axis (and 
therefore eventually fill in the O^-hole). Since we preserve the 
NO plume geometry, the total NO in the plume is proportional to 
its initial center line concentration and we use this value in 
this analysis. 

The results of these calculations are shown in Figure 25 
where the center line concentration of O3 at t = 150 sec are 

, ' p 

plotted as a function of the intensity of turbulence w 1 , and 

the initial center line concentration of NO. Recall that the 

initial center line concentration of O3 is always zero. For 

comparison with the previous calculation, also note that the 

initial center line concentration of NO assumed there is off-scale 

to the left in Figure 25. 


Examination of Figure 25 quickly shows that the penetration 
of 0^ to the plume axis increases with either an increase of the 
intensity of turbulence or a decrease of the amount of NO in the 


plume. However, 
10 -2 mysec 2 , the 
to about l/100th 
get significant 0 
period. The redu 
cent of present e 
would probably al 
operations in the 


2 -"5 

if we accept a value of w 1 between 10 J to 
total amount of NO in the plume must be reduced 
that used in the previous calculation in order to 
q penetration from the start of the post-vortex 
tion of NO production by jet engines to one per 
stimates would remove the Oo-hole anomaly, but it 
so resolve the problem of O3 depletion due to SST 
stratosphere. 1 


Viewed from another perspective, however, the results shown 
in Figure 25 begin to suggest a diurnal oscillation of 0^ penetra- 
tion of the exhaust plume. When the plume enters the dark sky and 
photodissociation of NO2 to produce 0 and NO ceases, the penetrat- 
ing O3 will transform the NO to NO2 thus reducing the NO amount 
(but not the total N0 X ). During the dark sky period we should 
expect Oo penetration of the plume, a calculation for which the 
model described here is suitable but which has not yet been made. 
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CONCLUSIONS AND RECOMMENDATIONS 


Given the objectives of this program to construct a coupled 
invariant diffusion and chemistry model and to exercise such a 
model sufficiently to show its applicability and advantages in 
real-world situations, this work has demonstrated the feasibil- 
ity of incorporating the stochastic nature of turbulent diffusion 
and chemistry in dynamic models, and has provided the first work- 
ing version of such a model. More importantly, perhaps, even 
this first and relatively simple version has revealed interactive 
effects between turbulent diffusion and chemical reactions which 
could not possibly be revealed by mean-value or first-order 
closure models. It seems clear to us that, on the one hand, the 
atmospheric chemists must reexamine their traditional assumptions 
of a well-mixed system in quasi -equilibrium, and on the other the 
atmospheric dynamicists must extend their considerations beyond 
the classical calculations of average values of pollutant concen- 
trations. The further exercise and. development of these second- 
order closure models can provide the necessary tools. 

We make this latter recommendation with a full awareness of 
the complexities of the chemistry of air pollution and the prolif- 
eration of the invariant model equations as multiple or chain 
reactions are introduced. However, this increasing complexity 
need not deter the development and use of these concepts, since 
they may first be used to analyze critical turbulent reactions, 
then to define areas where simpler models are quite adequate, and 
finally to provide simulation capabilities for those situations 
where first-order closure models are demonstrably inadequate. 

As a desirable prelude to further development of these models, 
such as their extension to coupled three-dimensional systems, the 
present model should be- subjected to a rigorous sensitivity analy- 
sis wherein the input variables of the turbulence field, the 
initial plume geometries, and the chemical reaction coefficients 
are systematically varied and the outputs of reaction and diffu- 
sion rates and the concentration distributions are tested for 
sensitivity to these input variations. Second, critical experi- 
ments, first in the laboratory and then in the atmosphere, should 
be designed and conducted to verify and validate not only the 
basic model output, such as mean values of concentration, but also 
the processes internal to the model's workings. In an atmospheric 
experiment, a minimum measurement program would require an array 
of towers oriented downwind from cross-wind line sources and 
equipped to measure the simultaneous means and fluctuations of at 
least two reacting chemical species (coming from well-defined 
sources) and the turbulent flux of these materials, all as a 
function of height. The conduct of such an experiment within 
the broader measurements program of the EPA/RAPS appears particu- 
larly desirable. 
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■ The conventional assumption of local uniform mixing of 
reactive chemical species is reexamined by derivation of the 
chemical reaction equations to include the effect of locally 
inhomogeneous mixtures on the reaction rates. Preliminary 
solutions of a simplified version of these equations show that 
inhomogeneities in reactant concentration generally tend to 
slow the reaction rate. Estimates of the relative roles of local 
diffusive mixing and chemical reactions in inhomogeneous 
mixtures show that there are several relatively fast photo- 
chemical reactions which may be limited by local diffusive 
mixing. In these cases, the reaction proceeds much more 
slowly than would be predicted if the reactants were uniformly 
mixed. 


I n developing either mathematical simulation models or 
laboratory chambers for the study of chemical reactions 
in the atmosphere, it has been generally assumed that the 
reacting materials are uniformly mixed. However, observations 
of the time history of concentrations of trace materials show 
quite clearly that uniformly mixed materials are the exception 
rather than tne rule in both air and water (Nickola et a!., 1970 
Singer et al., 1963; Csanady and Murthy, 1971). Local fluctua- 
tions of concentration are particularly significant during the 
early stages of atmospheric mixing, immediately following dis- 
charge of trace materials into the atmosphere, and when there 
are multiple point sources of pollutants. Our purpose here is to 
make a preliminary estimate of the importance of these fluctua- 
tions on atmospheric chemical reaction rates and determine, 
at least approximately, the relative roles of reaction rates and 
diffusive mixing in the control of atmospheric chemical re- 
actions. 


1 To whom correspondence should be addressed. 
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The importance of chemical reactions in the atmosphere 
has been increasingly recognized in the problems of air pollu- 
tion. These are probably most acute in dealing with photo- 
chemical smog formation (Worley, 1971). We have, therefore, 
drawn our examples from photochemistry, but we have not 
attempted to go beyond an examination of the possible im- 
portance of inhomogeneous mixing in these processes. 


Basic Chemistry Model 
We assume a bimolecular reaction 


a -{• —*■ y -h 8 (1) 

where a, /?, y , and 8 denote chemical species and that the re- 
action of a with /3 to form y and 8 is stoichiometric and is 
governed by equations of the form 


cM 

dt 

m 

dt 


-*[«][/?] 

-K[«m 


( 2 ) 


(/] denotes the molar concentration of the ith chemical 
species and K is the reaction rate constant in units of (sec 
mol/cm 3 ) -1 . It is convenient to transform the concentration 
terms in Equations 2 to dimensionless mass fractions, C ( , by 

PoCi = Mi[i] (3) 

where p 0 is the density of the mixture (g/cm 3 ) and is the 
molecular weight of the /th chemical species. Then the de- 
pletion rates for the « and j9 species may be written 

— = — K a C a C$ (4) 

dt 

and 

— s ~ ? C a C» (5) 

dr M a 


68 



APPENDIX A 


where K a - KpJM 0 and has dimensions (sec-ppm) -1 when C a 
and Cfi are expressed in parts per million (ppm) by wt. 

We may now examine the relative contributions of the means 
and fluctuations of C a and C 0 to the chemical reaction rate by 
assuming the time history of these quantities at a fixed point 
constitutes a stationary time series and that 


c a = C „ + ') 

Cp = Cp + Cp'j 


( 6 ) 


where the overbar indicates a time average and the prime 
indicates the instantaneous fluctuation about the average. 
Noting that CJ = Cp = Oand that 


ac< ^ 5>C< , be/ 

dt dt c )/ 


( 7 ) 


we obtain directly from Equations 4-6 

^ = -K a (C„Q + CVCV) (8) 


and 


dCp 

dr 


KM? 

M a 


(C a Cp + C a 'C/) 


( 9 ) 


where we have suppressed the dependence of K a on the tem- 
perature and pressure. (This analysis can be extended to in- 
clude the fluctuations of K a owing to significant fluctuations of 
temperature and pressure. For our present purposes, we shall 
assume an isothermal reaction at ambient pressure.) 

The role of concentration fluctuations in chemical reactions 
is immediately evident from Equations 8 or 9. The second- 
order correlation in the joint fluctuations of C a and C 0 either 
enhances the reaction rate (when the correlation is positive) or 
suppresses the reaction when C a 'Cp' is negative. Only when 
these fluctuations either do not exist or are uncorrelated is the 
average reaction rate governed by the average concentrations. 
As a simple example of the importance of this correlation 
term, imagine that the materials a and |3 pass the point of 
observation at different times — i.e., they are never in contact 
with each other. Values of C a and C 0 would be observed, but 
it is readily seen that C a Cp = — C a T C 0 ' in this case, a result 
which correctly predicts no chemical reaction. 

If we assume no diffusive mixing of the reacting materials, 
we may model the chemical reactions by noting that 


- — — l — — 

dt dt 


00 ) 


ac/ _ 

dC _ 

ac 

(ID 

dt 

dt ~ 

a / 

acycy „ , 
a/ * 

,dC 0 ' 

dt 

+ C ' &<v 
+ c * 

(12) 


and 


Performing the necessary operations and time-averaging, we 
get, repeating Equations 8 and 9, 


= -K a (c.c, +cycy) 


(13) 


^ = - ic ° c » + c <-' c e‘> no 

dt M a 


acyty = _ K I" 

dt “ lM a 


M*,- 


(Qcycy + c 0 c a '* + c a '*c p f ) + 


(CpCJCp’ + C a Cp 2 + c £ 




(15) 


7\C # « __ 

— - - -2 K a (C a C a 'Cp f + CpC a '* + Cp'C a '*) (16) 


= _ 2 ~^r- e c*c°’Ce' + c„cy» + c.'<y*) (i7) 

dt M a 

Equations 1 3-17 provide a closed set, except for the third-order 
correlation terms CJCp" 1 and C 0 'C a ' 2 . 

The appearance of the third-order correlations complicates 
the modeling problem very considerably since the statistical 
description now requires consideration of the distribution 
functions for C a and Cp. In an independent study, O’Brien 
(1971) has proceeded from Equations 13-17 by assuming the 
form of these distribution functions. Another approach, which 
we are pursuing, is to model the third-order correlations in 
terms of the second-order correlations. However, for our 
present purpose of determining whether or not the effects of 
inhomogeneous mixtures on chemical reactions may be 
significant, we may neglect the third-order correlations by 
assuming CJ and Cp are symmetrically distributed about C a 
and Cp, respectively. This assumption is, of course, untenable 
for more general cases but it does permit a solution of Equa- 
tions 13-17 by numerical techniques. 

For an initial test of the significance of inhomogeneities in 
chemical reactions, we assume a reaction box in which the 
initial concentration distributions of a and 0 are arbitrarily 
specified by Cy Cy CV 2 , cy 2 , and CJCp. As a further con- 
straint which isolates the chemical reaction process, we assume 
there is no mixing in the reaction vessel and no wall effects. 

As a reference case, let us assume a completely uniform 
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Figure 2. Chemical depletion of initially 
inhomogeneously mixed reactants (C'7C> 5 
=* 0.40) for various degrees of initial correla- 
tion between C a ' and Q>\ as measured by 

n _ CaCtf 
a3 



mixture of a and /3 — i.c., no fluctuations in concentration, 
and initially C a — Cp — C 0 . The predicted values 
of C a are shown in Figure 1 as a function of time normalized 
by the reaction rate constant. As can be seen, the reaction 
proceeds to exhaustion of the reacting materials. 

Now let us assume that a and /3 are initially inhomogene- 
ously mixed but that there is no initial correlation between 
C a r and cy--i.e., initially C a 'Cf 0. Asa measure of these 
fluctuations, we take C t ' 2 /Q 2 = 0.2, 0.4, 0.6, 0.8, and 1.0. 
The results of these calculations are also shown in Figure 1, 
and it is immediately evident that any inhomogeneities 
operate to suppress the chemical reaction rate and to stop it 
completely before the reacting materials are exhausted. 
Mathematically, the model predicts that, in the absence of mix- 
ing, initial inhomogeneities operate to produce values of 
C a 'Cf which eventually become equal to — C a Cp and the re- 
action ceases. Physically, the local reactions have everywhere 
proceeded to exhaustion of one of the reactants, leaving a 
residue of the other reactant and products at that site. 

It is of special interest to note, from Equation 15, that the 
suppression of the reaction rate by C a 'Cy depends only on 
one of the reactants being nonuniformly distributed initially. 
A negative rate of change of C a 'Cy can be generated by non- 
zer o va lues of either C a ' 2 or Cy 2 , since the terms C^C a ' 2 and 
C a Cy 2 are positive definite. The presence of concentration 
inhomogeneities in one of the reactants generates inhomogene- 
ities in the other. 

The effect of an initial correlation between C a ' and Cy may 
now be examined by assigning initial nonzero values to C a 'Cy , 
C a /2 , and Cy 2 . To illustrate this effect, we have chosen 
C7VC, 5 = C/yCV = 0.4 and CVCZ/CCV’CT’)" 2 - 
+ 1.0, +0.5, 0.0, -0.5, and -1.0 where C„'C fl '/(C a ,5 C s ' 5 )'' 2 
== Ra.ti, the ordinary correlation coefficient. The resulting 
predictions of C(C 0 are shown in Figure 2 and are again com- 
pared with the uniformly mixed case. As might have been ex- 
pected, initial positive correlation accelerated the reaction 
rate, but only when this initial positive correlation was per- 
fect did the reaction go to exhaustion of the reacting materials. 
In this case, although there were concentration fluctuations, 
stoichiometrically equal amounts of a and /3 were initially 
placed in each local volume. In all other cases, the reaction 
was again halted when one of the reactants was exhausted 
locally, leaving a residue of the other reactants and products 
of the reaction. 

The combined effects of initial inhomogeneities and correla- 
tions between the fluctuations are summarized in Figure 3 by 
plotting the depletion of C during the first normalized time 
step as a function of R a i<3 and C' 2 /C 0 2 . The effect of the mag- 



Figure 3. Joint effect of initial correlation and inhomogeneity on de- 
pletion of reacting materials at Kt = 1.0 


nitude of the fluctuations, as measured by C'\ reverses as one 
goes from large positive toward small positive and negative 
values of C a 'Cy . 

These results point toward an important role for fluctua- 
tions of concentration in controlling chemical reaction rates. 
For example, if two reacting materials are discharged si- 
multaneously from a point source, during their initial mixing 
with the atmosphere their concentration fluctuations should 
be large and positively correlated. We would then expect, on 
the basis of this effect, that their reaction rate would be con- 
siderably faster than if they were uniformly mixed from the 
start. The emission of hydrocarbon and NO* from auto ex- 
hausts is a case in point. Discharge of S0 2 and particulate 
matter from power plant stacks is another. 

On the other hand, if reacting materials are randomly 
mixed or if positive fluctuations in one are associated with 
negative fluctuations in the other, the reaction should be 
suppressed, compared to the uniformly mixed case. Both of 
these cases could be important, but their true importance de- 
pends critically on the rate at which atmospheric diffusion 
tends to mix chemical species and, hence, to diminish these 
fluctuations, as compared with the rate of chemical reaction 
produced by the concentration fluctuations. 

Estimates of Local Mixing Rates in the Atmosphere 

The only way in which the correlation C a 'Cy can be elim- 
inated, if it exists, in given flow situations is by the process 
of molecular diffusion. To estimate the rate at which this can 
occur, we may write the expressions for the contribution of 
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molecular diffusion to the time rate of change of C a ' and Cp'. 
They are 


W'CJ = D a 


( d - c A -D. 

\ d/ / diff 

(*£A 

\ <>/ /dirr by&yt 


&c a f 

dy^yi 


(18) 


(19) 


Multiplying Equations 18 and 19 by Cp and C a \ respectively, 
adding, and time-averaging gives 


= D 

\ df /diff 


VCJCfi' 2 D dc a ' (SCp 


( 20 ) 


(We have assumed D a ^ Dp — D, consistent with the assump- 
tion M a Mp. See O'Brien (1971) for a discussion of this 
assumption.) The first term on the right-hand side of Equation 
20 is nondissipative— i.e., it measures the transfer of the CJCp 
correlation by gradients in the value of this correlation within 
the field. The second term is dissipative — i.e., it measures the 
local diminution of C a 'Cp by the action of molecular diffu- 
sion. The appropriate expression for this term is 


2 D <>£V = 2 D C a 'Cp' 
dy dy X 2 


( 21 ) 


In this expression, the dissipative scale length X must be chosen 
as it is chosen for the calculation of other turbulent correla- 
tions when performing calculations of the structure of turbu- 
lence (Donaldson, 1969). 

For such calculations, 1/X 2 is given approximately by 

1 ^ 0.05 p„g 
X 2 ~ p 0 A 

where p 0 is the atmospheric density, q 2 = u" 1 -F v ' 2 + h>' 2 , p 0 
is the molecular viscosity of air, and A is a length scale related 
to the integral scale of the atmosphere, and is of the order of 
1000 cm in the earth’s boundary layer. 

If we choose typical values of the parameters involved in 
evaluating the magnitude of the expression for ( dC„'Cy/d/)diff 
given in Equation 20, we have 


A = 1000 cm 
Po = 10- 3 g/cm 3 
p 0 - 1.7 X 10“ 4 g/cm-sec 
D = 1.7 X 10~ 1 cm 2 /sec 
q = 30 cm/sec 

These numbers give for the magnitude of the dissipative scale 
X = 10 cm 


From this result, we obtain, finally, 

2 D 

= 3.4 x 10“ a sec -1 

X 2 

The rate of destruction of CjCp by molecular diffusion is of 
the order of magnitude of 3.4 X 10* 3 C a 'Cy and is tending to 
drive C a ’Cp to zero. If this diffusion dissipation of C a ’Cp is 
dominant, CJCp will remain close to zero and the chemical 
reaction will proceed according to the product of the mean con- 
centrations. In this case, the chemical reaction rate is controlled 
by the reaction rate constant and the mean concentrations in 
the traditional way. 

On the other hand, if the change of CJCp is dominated by 
chemical reactions — i.e., the reaction proceeds more rapidly 
than species diffusion, CJCp' will tend to the value — C a C$ 


Table I. Reaction Rate Constant K a for Various Photochemical 
Reactions and Associated Estimates of N Using C a = 1 ppm 


Reaction 


K (ppm-sec)~ 


I. 

0 3 + NO = N0 2 + 0 2 

8.3 X 10-* 

4 

2. 

N0 2 ~h 0 3 — NOj T O a 

1.7 X 10“ 5 

2 X 10 2 

3. 

N0 3 -F NO = 2N0 2 

4.8 

7 X 10“< 

4. 

NO + HO, = NO, -F OH 

1 .7 X 10- 1 

2 X 10~ 2 

5. 

OH -F 0 3 = HO, -F 0 2 

1.7 

2 X 10“ 3 

6. 

OH -F CO = H -F C0 2 

5.0 X 10“ 2 

7 X 10- 3 

7. 

ch 3 o 2 -f no = ch 3 o -f no 2 

1.7 

2 X 10“ 3 

8. 

c 2 h 3 o 2 -f no = 




c 2 h 3 o -f no 2 

1.7 

2 X 10“ 3 

9. 

c 2 h<o 2 -f no = 




CHaCHO -F NO, 

1.7 

2 X 10“ 3 

10. 

CH a O -F 0 2 = HCHO + HO, 

1.7 

2 X 10“ 3 

11. 

C 3 H fi -F O = CH 3 -F C 2 H 3 0 

6.0 X 10" 1 

6 X 10“ 3 

12 . 

c 3 h 6 -f o 3 = 




HCHO -F C 2 H,0 2 

8.3 X 10" 3 

4 X 10" 1 

13. 

C 3 He -F 0 2 — CH a O “F C 2 H 3 0 

1.7 X 10“ 2 

2 X 10" 1 

14. 

c 3 h 6 + ho 2 = 




ch 3 o -f ch 3 cho 

3.4 X 10~ 2 

lO” 1 

15. 

c 2 h 3 o -f m = 




CH 3 + CO + M 

1.7 X 10“» 

2 X 10“ 2 


and the reaction rate will be suppressed. In this case, the re- 
action is controlled by the rate of species mixing and will de- 
pend on parameters other than K a and C a Cp. 

We may estimate which two-body reactions will proceed as 
though C'aCp 0— i.e., in the usual manner, and which will 
be modified by having values of |C a 'Cy| of the same order as 
|C a Q| by forming the ratio 



When N 1.0, CJCp will tend to zero and the reaction will 
be controlled by the reaction rate constant and the mean con- 
centrations; when TV <$C TO, CJCp will tend to and remain 
close to -C a C & and the reaction will proceed at a rate deter- 
mined largely by the rate at which one reactant can be mixed 
with another and will depend on the scale of the patches of 
unmixed reactants. 

Typical values of K a for various reactions which enter into 
the photochemical chains are listed in Table I along with esti- 
mates of N. For these two-body reactions, only the first and 
second are sufficiently slow for conventional kinetic models 
to apply. The propylene reactions with 0 3 , 0 2 , and H0 2 (num- 
bers 12, 13, and 14 in Table I) tend to represent a transition 
stage between diffusive mixing control and chemical reaction 
control of the reaction rate. The remaining reactions are all 
clearly diffusion-limited in inhomogeneous mixtures and 
should proceed at a rate which is much slower than conven- 
tional chemical kinetics would predict. 

Conclusion 

These results indicate there are clearly reactions in the photo- 
chemical chain which will be suppressed by the inability of the 
atmosphere to mix the reacting materials rapidly enough to 
prevent serious local depletion of one of the reacting materials. 
In these cases, conventional models of the reaction will tend to 
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seriously overestimate the reaction rate, and, therefore, the 
production rate of the chemical species which enters into the 
next reaction in the photochemical chain. On the other hand, 
the enhancement of the reaction rate for two materials 
emanating from a common source and, therefore, occupying 
the same volume of the atmosphere, during the initial period 
of incomplete mixing, may also represent a significant de- 
parture from conventional simulation models. 

It is hoped that this brief and necessarily incomplete dis- 
cussion will serve to demonstrate the importance of turbulent 
fluctuations of concentrations in atmospheric chemical re- 
actions. Consideration of these effects in refining simulation 
models of these reactions appears to be important. 

Nomenclature 

Ci = mass fraction of ith chemical species, ppm 
D t = molecular diffusion coefficient for ith chemical 
species, cm 2 /sec 

K = chemical reaction rate constant, cm 3 /sec-mol 
/f a , Kp = chemical reaction rate constants, 1/ppm-sec 
Mi — molecular weight of ith chemical species, g/mol 
N = nondimensional ratio of characteristic times 
q 2 = u' 2 + o' 1 -j- w' 2 , cm 2 /sec 2 

R a p = ordinary second-order correlation coefficient 
/ = time, sec 

u’ , o', w' = orthogonal components of turbulent fluid 
motion, cm/sec 

— length along ith direction of a cartesian coordinate 
system, cm 


{ ], “ = averaged quantity 

' = departure from the average of the primed quantity 

Greek Letters 

a,f},y,8 = chemical species 

\ — dissipation scale length, cm 

A = macroscale of atmospheric turbulence, cm 
n 0 = dynamic viscosity for air, g/cm-sec 
p 0 ~ fluid density, g/cm 3 
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CHEMICAL REACTIONS IN INHOMOGENEOUS MIXTURES: THE EFFECT OF THE 

SCALE OF TURBULENT MIXING 


Coleman duP. Donaldson* and Glenn R. Hilst** 


ABSTRACT 

Recent studies by O’Brien [1] and the authors of this paper [2] 
have provided a theoretical framework for the assessment of 
chemical reaction rates when the reactants are embedded in a 
turbulent fluid and are inhomogeneously mixed. The results of 
these studies, which are reviewed here, point towards a profound 
effect on chemical production and depletion rates when the charac- 
teristic reaction time, as measured by the product of the chemical 
kinetic reaction rate constants and the average and fluctuating 
concentrations of the reactions, is short compared with the 
characteristic molecular diffusion time. The latter time is 
measured by the ratio of the molecular diffusion coefficient and 
the square of the dissipation scale length, and is, therefore, 
dependent upon the scale of the turbulent motions. Both the fact 
of inhomogeneous mixtures and this dependence upon turbulent 
scales of motion pose significant problems when extending labora- 
tory results to other scales of motion, such as the free atmosphere. 

These theoretical results, which are partially substantiated 
by observations, point towards the need for simultaneous measure- 
ments of turbulence and chemical reaction rates over a range of 
turbulence scales and reaction rate constants. If substantiated 
by such new experimental measurements, the theoretical results 
point towards a clear requirement for joint- consideration of the 
chemical reactions and the scale of turbulence in such diverse but 
critical problems as the design of large combustion apparatus and 
the calculation of photochemical reactions in the atmosphere. 


INTRODUCTION 

Although the effects of inhomogeneous mixing of reacting 
chemical species on the reaction rate, as measured by either the 
depletion of reacting species or the production of new species, have 
been recognized for at least ten years [ 3 ], methods for accounting 
for this effect have only recently emerged [1,2]. Neither of 
these methods are as yet fully developed, but they are sufficiently 
advanced that we may make some preliminary estimates of the 
situations under which the effects of inhomogeneous mixing will be 
pronounced or perhaps even completely dominate the reaction. 

In particular, we find for .the case of an irreversible two- 
body reaction at constant temperature that the following limita- 
tions are imposed by inhomogeneous mixing of either or both of the 
reacting species: 
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• If the chemical reaction rate is slow compared with the 
molecular diffusion rate, no effect is noticed, and the 
reaction proceeds according to conventional chemical kinetics . 

• If the chemical reaction rate is fast compared with the 
molecular diffusion rate, the reaction rate is limited by 
the diffusive mixing rate, tending, in the limit of very 
slow diffusive mixing, to zero before the reactants are 
exhausted . 

A large number of reactions in . combustion processes and photo- 
chemical smog formation fall within this latter category. It is, 
therefore, of considerable interest to investigate further just 
how much the reaction rate is curtailed by inhomogeneous mixing 
under such circumstances. In the following pages, we derive the 
basic equations for prediction of the joint effects of chemical 
reactions and molecular diffusion, examine the effects of the 
dissipation scale length of the turbulent motions, and identify, 
on a preliminary basis, the two-body reactions inherent in photo- 
chemical smog formation for which inhomogeneous mixing is a 
limiting condition. 


MODELING OF CHEMICALLY REACTING FLOWS 


For most computations of chemically reacting turbulent flows, ' 
it has been customary for engineers to proceed with the calculation 
according to the following scheme. First, the engineer develops by 
some method (mixing length, eddy diffusivity, or other method) 
equations for the time-averaged or mean values of the concentrations 
of the reacting species of interest (say, species a and P ) at 
each point in the turbulent flow under consideration. He also 
obtains an equation for the mean value' of the temperature that is 
expected at each point in this flow. It is then customary, if the 
equations that generally govern the reaction between a and P , are 



-k l C cx C P 


( 1 ) 


DCp 

Dt 


-k 0 C C R 
2 a p 


( 2 ) 


to assume that valid equations for the time rates of change of the 
mean values of the mass fractions of a and P are 

( 3 ) 

(4) 

In these equations, C^ and C^ are the time-averaged mass 
fractions of the two species and k^ and k^ are the reaction 
rates k^ and k^ evaluated at the mean temperature T , i.e.. 


DC a 

= -k, C C„ 
Dt lap 


dc 6 

Dt~ = ~ k 2 C a C p 
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k 1 = k 1 (T) and k ? = k 2 (T) . 

Although equations such as (3) and (4) are used extensively at 
the present time, it is not difficult to show that they are 
incorrect when reaction rates are fast and the scale of the turbu- 
lence is large. This may be done by considering the proper forms 
of Eqs . (1) and (2) when they are averaged. The well-known results 
are* 


DC 

a 

Dt 


f 




~ 


- k i( 

C C Q + 
v a p 

C'C') - 
a p/ 

C k'Cl + 
alp 

C p k i C a + 

Vc i C } C ’ 
K l° a L p 

(5) 

DC r 
_E = 
Dt 


< _ 


r- 


-| 


-k 2 ( 

C C a + 
, a P 

C'C') - 
a $) 

C k »C» + 
a 2 p 

% k 2 C a + 

k * c } r 1 

K 2 a°p 

(6) 


To demonstrate the character of these equations, let us discuss 
them under the assumption that k-^ = k£ = 0 . Equations (5) and 
(6) then reduce to 


and 


DC / 

Dt 5 ’ - k l( C a C p + C a C p) 


DC 


BiT ' - k 2 ( C a C p + °i C fs) 


— / 


(7) 


( 8 ) 


It is clear from these equations that, if one wishes to calculate 
the reaction of a with (3 , it will be nec essa ry to have an 
equation for the second-order correlation C ? CA unless one can 

_ __ CL p 

show that f° r P ar ^i cu ^ ar flow in question. The 

conditions required for C ' C ' « C can be derived in the follow- 
er p up 

ing way. First, by following the method used by Reynolds for the 
derivation of the equation for the turbulent stress tensor, one 
finds the fol lowi ng equation for the substantive derivative of the 

correlation C'C' :** 
a P 


DC a C p 

Dt 


DC 


Dt 


<£k) 


1 chem 


- u J C 


a p, j 


. - u J r C ' C . - u J 'C'C 


P 


+ % (C^) >mn - 2 & g 11 ( C 


a p ,J 


a,m p,n 


(9) 


* For a discussion of these equations that is related to the 

present treatment, reference should be made to O'Brien [1] which 
was published after this work on the modeling of chemically 
reacting turbulent flows was started. 

** For the purposes of this illustrative discussion, the flow is 
treated as incompressible. 

T The notation is that of general tensor analyses. g mn is the 
contravariant form of the metric tensor g mn . 
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where the term (DC’C’/Dt) , is the contribution of chemical 
a p chem 

kinetics alone to the substantive derivative of C*C' . This 

ex p 

expression can be found from Eqs . (1) and (2), and is 


PC q C f3 

1 = -kj 

f ChCl 2 + C Q C 'C ' 

+ C'C' 2 ) 

Dt 

/chem 1 

( a £ p a P 

a p / 


-k~| 

[ C C'C' + C D C' 2 

+ C' 2 C’ ) 


2' 

^ a a p pa 

a p / 


( 10 ) 


It is instructive to discuss the behavior of the correlation C MI ’ 
for the case of turbulent reactions in the absence of any p 
appreciably large gradients. In this case, Eq. (9) becomes 


DC a C (3 _ 
Dt 


DC a C p 

Dt 


chem 


- 2c0g mn ( C 


' C ' ) 

a,m p,n' 


( 11 ) 


The second te rm on the right of Eq . (11) is the destruction of the 
correlation C^CL by the action of molecular diffusion. In line 
with our previous work [5] , we will model this term by means of a 
diffusion scale length X so that Eq . (11) becomes 


DC f C ’ 
a ft = 

Dt 


dc ; c m 


Dt 


- 2 <& 


C 1 Cq 
a ft 


/chem 


( 12 ) 


The diffusion term in this equation is such that 
approach zero with a characteristic time that is 


T dif f 


\ 2 /2& 


C'C' tends to 
a ft 

(13) 


What is the overall effect of the first term on the right-hand side 
of Eq . (12)? The effect is difficult to see from an inspection of 
Eq . (11), but we may derive an expression for what this term _ 

accomplishes from Eqs. (7) and (8). First, multiply (7) by Co and 
(8) by C a and then add the resulting equations. The result is 

( DC C A \ . _ w 

' '( kA + * °«°b) <14) 


This equation can be interpreted by saying that the eff ect of c hem- 
istry alone is to drive C Q Cp to the negative of C^C^ (or C^CI 
to the negative of C a Cp) with a characteristic time p p 


T 


chem 


k l C p + k 2 C a 


(15) 


Equation (14) states that the reaction between a and (3 will 
always stop, i.e., c a c p + will become zero, short of the 

exhaustion of a or ft unless a and ft are perfectly mixed wherever 
they occur in the turbulent flow under consideration. The physical 
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reason for this is that, in the absence of diffusion, if a and P 
are not perfectly mixed to start with, the final state of the gas 
in any volume element will be a and products, P and products, 
a alone, or P alone, but never any region containing both a and 
P . It is easy to see that, no matter what the values taken on by 
and Cp are as a function of time, if C Q is never nonzero 

when Co is nonzero and vice versa so t hat n o reaction is possible, 
it is mathematically true that CC Q + C'C' = 0 . Thus, Eqs . (7) 

U p U p 

and (8) state that no reactions are possible as required by the 
physics of the problem. 


An actual example may make the meaning of C^C^ more clear. 

Consider that the flow of material by a given point is such that 
alternate blobs of a and p pass the point. Let us suppose that 
half the time the flow is all a and half the time it is all p . 
The resulting concentrations are sketched in Figure 1. If this 
pattern keeps repeating, the average values of C a and are 

obviously C Q = 1/2 and Cg= 1/2 . Whenever the flow is all cl y 

=+1/2 and C£ =-1/2. Whenever the flow is all P , =-1/2 

and C' = +1/2 . We find then that the average value of C 1 C' 

P _ _ CL p 

must be C'C 1 = -l/J| . Since C T C* = C , no reaction is possible 
a p a p a p K 

according to Eqs. (7) and (8) and obviously no reaction should 
occur . 


THE EFFECT OF SCALE LENGTH 


We may now return to Eq. (12). If, in this equation, the 
scale A is small enough and the reaction rates are slow enough, 
the second term on the right-hand side of the equation will be 

dominant and the flow will be such that C'Ci is always almost 

a p 

zero* This means that molecular diffusion is always fast enough 
to keep the two species well mixed. On the other hand, if the 
reaction rates are very fast and X is very large, the fi rst term 
on the right-hand side of Eq . (12) will be dominant and C^C^ will 

tend to be approximately equal to -C^Cp and the two species will be 

poorly mixed. The rate of removal from the flow of a and p by 
reaction will then not be governed by reaction rates but will be 
limited by molecular diffusion. To put these notions into quanti- 
tative form, let us consider the ratio of the two characteristic 
times 


N 


T dif f 
T chem 



k l C (3 


k 2 c 


a) 


( 16 ) 


and a contact index 


I 


C C + C'C' C'C' 

a P a p _ 1 + a P 


C C R 
a P 


c C A 

a P 


( 17 ) 


We note that if N is much smaller than one, diffusion will be 
very rapid and the two species a and p will be in intimate 


78 



APPENDIX B 


258 Donaldson and Hilst: Chemical Reactions and Turbulent Mixing 


contact with each other. In this case C!C'/C C_ will be small 

a p up 

and the contact index will approach one. If, on the othe r hand, 
N is much larger than one, mixing will be poor and C 1 C ’ will 

Q p 

approach -C Cg . The contact index will then approach zero. 

In this case xhe reaction will be diffusion-limited. 


In many laboratory flows, the dissipative or diffusive scale 
of turbulenc e is very small and N is, indeed, small so that the 
neglect of in the kinetic equations is permissible. On the 

other hand, if the laboratory experiment is just increased in size, 
holding all other parameters such as velocity, temperature, etc., 
constant, one soon finds that the character of the flow changes. 
This may be seen by examining the expression for the dimensionless 
quantity N in more detail. 


Let us assume the diffusive scale of a turbulent flow is of 
the order of the dissipative scale so that we may relate X to the 
integral scale length of the turbulence A 1 by (cf. Ref. [5]) 

x 2 = A 2 /(a + bpqAj/n.) (18) 

where a and b are constants and pqA_,/p is the turbulent 
Reynolds number. Substituting this expression into Eq . (16) gives 


N 


n i/ k i c e 

2£>\a + 


+ k 2 C a \ 
bpqA 1 /p.j 


(19) 


For relatively high Reynolds numbers, this 

N = • - 75 : * — (k-iCft 

2 b q \ 1 P 

If an experiment is performed in the laboratory and a value of N 
for this experiment is determined or estimated and is found to be 
small compared to one, then we know that the diffusive mixing of 
the flow is such that the species a and (3 are in contact. The 
reaction rate of these species is then chemically controlled. Now 
if the apparatus is just scaled up in size, all other things being 
equal, N will increase linearly with size since the scale 
increases linearly with the size of the apparatus. When the scale 
has been increased sufficiently, so that N is no longer very 
small compared to one, the nature of the flow in the device must 
change, for the species a and (3 will no longer be in intimate 
contact at equivalent positions in the apparatus. 


expression becomes 
* k 2 f a) 


The turbulent atmospheric boundary layer is a good example of 
a fl ow in which it is essential to keep track of the correlation 

if one is to be able to make sense of the reaction of species 
which are introduced into the flow. To demonstrate this, we list 
in the table some of the second-order reactions responsible for the 
production of photochemical smog. We have also listed in this 
table the reaction rate recommended for each reaction [4] and an 
estimate of the number N for each reaction if it occurs in the 
atmospheric boundary layer where a typical value for X is 10 
centimeters. It is interesting to note that it is, in general. 
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those reactions listed in the table for which N is greater than 
one that investigators have found to proceed more slowly than pre- 
dicted by formulas such as Eqs . (3) and (4) when the reaction rate 
determined from laboratory experiments is used. This difficulty 
has led some investigators to search for other chemical reactions 
that might be considered which would explain this discrepancy. 

CONCLUSION 

It certainly appears unwise to follow this course until such 
time as one has at least developed a viable scheme for properly 
computing turbulent reacting flows. It is the authors* opinion 
that an acceptable method of computing such flows can be developed 
through the use of second-order correlation equations such as 
Eqs. (9) and (10). Methods of modeling the third-order correlations 
that appear in these equations can be found that are similar to 
those used to study the generation of turbulence and turbulent 
transport [5]- The development of a viable method for computing 
chemically reacting turbulent flows according to such a scheme is 
under active development by the authors. It is important to note 
in this connection that it is essential in developing this general 
mbthod to consider fluctuations in density and in the reaction rate 
constants when the chemical rate equations are considered. 


a , b 


CJ 

l 


Jdr 

I 


q 

A 

\ 

M* 

P 

T 


NOMENCLATURE 

= constants 

= concentration of subscript species, expressed as a 
mass fraction 

= molecular diffusion coefficient 
= contact index (Eq. (17)) 

= chemical reaction rate constants 

= dimensionless ratio of characteristic times for mole- 
cular diffusion and chemical reaction 

= rms value of turbulent kinetic energy 

= integral scale length of turbulence 

= dissipative scale length 

= viscosity 

= density of the fluid 
= characteristic time 
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Some Second-Order Reactions Responsible for Photochemical Smog [4] 


Reaction 

k (ppm-sec) - * 1 

N 

0 3 + NO = N0 2 + 0 2 * 

8.3 * 10 -4 

0.25 * 

no 2 + O 3 = no 3 + o 2 

1.7 x 10 -5 

5.0 x 10~ 3 

N0 3 + NO = 2N0 2 

4.8 

1.4 x 10 3 

NO + H0 2 = N0 2 + OH 

1.7 x 10' 1 

50.0 

OH + 0 3 = H0 2 + 0 2 

1.7 

5.0 x 10 2 

OH + CO = H + C0 2 

5.0 x 10 -2 

1.5 x 10 2 

ch 3 o 2 + no = ch 3 o + no 2 

1.7 

5.0 x 10 2 

C 2 H 3 0 2 + NO = C 2 H 3 0 + N0 2 

1.7 

5.0 x 10 2 

C 2 H 4 0 2 + NO = CH^CKO + N0 2 

1.7 

5.0 x 10 2 

CH 3 0 + 0 2 = HCHO + H0 2 

1.7 

5.0 x 10 2 

C + 0 = CH 0 + C H o 0 
3 b 3 2 3 

6.0 x 10" 1 

1.8 x 10 2 

C 3 Hg + 0 3 = HCHO + C 2 H 4 0 2 

8.3 x 10" 3 

2.5 

C^H. + 0 o = CH o 0 + C ~H o 0 
3 d 2 3 2 3 

1.7 x 10 -2 

5.0 

C -.Hs + H0 o = CH o 0 + CH.CHO 
3 d 2 3 3 

3.4 x 10 -2 

10.0 

c 2 h 3 o + m = ch 3 + co + m 

1.7 x 10 _1 

50.0 

#Those reactions for which 

N is small compared to 

one are those 

which can be treated using 

mean quantities in the 

basic equations 


of chemical change, i.e., correlations in fluctuating quantities 
may be neglected. 
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C 



Figure 1. Simple problem illustrating C'Ci = -C C when no 
reactions are possible op a p 
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